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I.  Introduction 


In  1993,  the  US  Department  of  Defense  announced  the  FY94  competition  for  the 
Augmentation  Awards  for  Science  and  Engineering  Research  Training  program.  In  order  to  be 
eligible  for  one  of  these  awards,  a  faculty  investigator  must  have  had  an  existing  DOD-funded 
research  grant  or  contract.  In  the  case  of  this  particular  award,  the  linked  DOD-funded  research 
grant  was  USAF  Grant  No.  F49620-93-1-0270  entitled  “Mechanics  of  Wide-Spread  Fatigue 
Damage  and  Life-Enhancement  Methodologies  for  Aging  Aerospace  Structures.”  Both  this  grant 
and  the  linked  grant  are  under  sponsorship  of  the  Air  Force  Office  of  Scientific  Research. 

The  distinction  between  this  grant  and  the  linked  grant  is  that  funds  provided  under  this 
grant  could  be  used  only  for  salaries  for  graduate  and  undergraduate  students,  tuition/fees  for 
graduate  students  only,  and  incidental  costs  associated  with  the  activities  of  the  students 
supported  under  this  grant.  Students  supported  under  this  grant  must  be  U.S.  citizens  or  nationals 
of  the  United  States  whose  research  activity  was  directly  related  to  the  linked  grant,  and  who  are 
in  addition  to  students  who  were  normally  accommodated  by  the  linked  grant. 


II.  Students  Supported  under  this  Grant 

During  the  course  of  this  grant,  three  students  were  supported  from  funds  provided  by  it. 
These  students  were: 

a.  Cummings.  Robert:  period  of  support:  1/3/95-6/9/95.  Mr.  Cummings  was  a  Ph.D. 
student  in  Materials  Engineering.  At  the  time  of  his  entry  into  that  program,  and  the  beginning  of 
support  under  this  grant,  Mr.  Cummings  was  employed  as  Materials  Engineer,  Warner  Robins  Air 
Logistics  Center,  Georgia.  Mr.  Cummings  did  not  prepare  any  technical  publications  while 
supported. 

b.  Jansen.  Randv:  period  of  support:  3/28/94-10/27/94.  Mr.  Jensen  was  a  Ph.D.  student  in 
Aerospace  Engineering.  At  the  time  of  his  entry  into  that  program,  and  the  beginning  of  his 
support  under  this  grant,  Mr.  Jensen  was  employed  as  a  structural  engineer  at  the  Warner  Robins 
Air  Logistics  Center,  Georgia.  Mr  Jansen  participated  in  the  research  described  in  the  enclosed 
report,  “Composite  Patch  Repairs  of  Cracked  Metallic  Structures,”  and  the  Progress  Report, 
“Mechanics  of  Widespread  Fatigue  Damage  and  Life-Enhancement  Methodologies  for  Aging 
Aircraft  Structures,”  which  was  submitted  in  connection  with  grant  no.  F49620-93- 1-0270, 
described  above,  also  enclosed. 


c.  Visser,  Bradley:  period  of  support:  6/21/93-12/9/94.  At  the  time  Mr.  Visser  received 
initial  support  under  this  grant,  he  was  a  graduate  student  (MS  level)  in  the  School  of  Aerospace 
Engineering  at  Georgia  Tech.  He  received  his  Master  of  Science  degree  in  June  1994  and 
entered  the  Ph.D.  program  in  Aerospace  Engineering.  While  employed  under  this  grant,  Mr. 
Visser  prepared  the  enclosed  report  “Summary  of  Research:  Alternating  Method  for  Bending  of 
Plates.”  Mr.  Visser  also  participated  in  the  research  described  in  the  enclosed  report, 
“Composite  Patch  Repairs  of  Cracked  Metallic  Structures,”  and  the  Progress  Report, 
“Mechanics  of  Widespread  Fatigue  Damage  and  Life-Enhancement  Methodologies  for  Aging 
Aircraft  Structures,”  which  was  submitted  in  connection  with  grant  no.  F49620-93- 1-0270, 
described  above,  also  enclosed. 

III.  Report:  “Composite  Patch  Repairs  of  Cracked  Metallic  Structures” 

This  report  appears  on  the  following  pages  numbered  1-57.  " 

IV.  Progress  Report  -  “Mechanics  of  Widespread  Fatigue  Damage  and  Life  Enhancement 
Methodologies  for  Aging  Aerospace  Structures.” 

This  two-page  report  follows  Part  III. 

V.  Report  -  “Summary  of  Research  -  Alternating  Method  for  Bending  of  Plates.” 


This  section  is  last. 
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CHAPTER  III 

COMPOSITE  PATCH  REPAIRS  OF  CRACKED  METALLIC  STRUCTURES 

3.1  Introduction 

In  order  to  enhance  the  life  of  cracked  structures,  the  repair  technique  which  uses  adhesively  bonded  com¬ 
posite  patches  is  being  considered  as  a  cost-effective  and  reliable  method[Baker  and  Jones  (1988)].  To  study 
the  effect  of  the  design  parameters  (including  the  size,  thickness,  and  material  properties  of  the  composite 
patch,  and  the  mechanical  properties  and  thickness  of  the  adhesive  layer)  on  the  crack-tip  stress  intensity 
factors,  finite  element  alternating  method  has  been  developed  for  the  studies  of  composite  patching  prob¬ 
lems.  The  section  summaries  the  procedure  used  for  the  analyses  of  composite  patching  problems  with 
numerical  examples  which  illustrate  the  effect  of  patching. 

The  traction  in  the  adhesive  layer  must  be  solved  for  first  in  order  to  solve  the  patching  problem.  Con¬ 
sider  the  following  two  problems.  The  first  problem  is  the  cracked  structure.  The  patch  is  replaced  by  the 
undermined  adhesive  traction  that  are  applied  on  the  area  where  the  patch  is  applied.  The  second  problem 
is  the  composite  patch  itself  with  the  adhesive  traction  applied  on  it.  Together  with  the  finite  element  alter¬ 
nating  method  for  the  cracked  structures,  two  different  approaches  had  been  used  to  solve  for  the  adhesive 
traction,  one  implicit  method  [see  Park,  Ogiso,  and  Atluri  (1992)]  and  one  explicit  method  [see  Nagaswamy, 
Pipkins,  and  Atluri  (1995)].  In  the  implicit  method,  the  stiffness  of  the  cracked  structure,  at  the  given  exter¬ 
nal  loading,  with  respect  to  the  traction  loading  is  determined  using  the  finite  element  alternating  method. 
A  linear  equation  system  can  be  obtained  for  the  adhesive  traction  by  using  the  compatibility  condition  be¬ 
tween  the  sheet  and  the  patch.  After  solve  for  the  adhesive  traction,  the  stress  intensity  factors  can  be  easily 
determined  through  the  analytical  solutions  used  in  the  finite  element  alternating  methods  for  the  cracked 
structure.  More  details  about  the  implicit  method  can  be  found  in  Park,  Ogiso,  and  Atluri  (1992).  In  the 
alternative  explicit  method,  the  adhesive  traction  is  solved  by  solving  the  patched  structure  using  an  ordinary 
finite  element  method,  with  the  crack  being  modeled  explicitly.  However,  the  crack  tip  is  not  modeled  in 
detail  since  the  objective  is  only  to  obtain  the  adhesive  traction.  Once  the  adhesive  traction  is  obtained, 
the  stress  intensity  factors  can  be  determined  accurately  by  analyzing  the  unpatched  cracked  structure  with 
adhesive  traction  applied  at  the  location  of  the  patch,  using  the  finite  element  alternating  method. 

3.  LI  Analysis  of  repaired  cracks 

In  this  section,  the  problems  of  composite-patch  repair  of  (i)  center  and  edge-cracked  panels  loaded  in  the 
far  field;  and  (ii)  cracks  emanating  from  pin-loaded  fastener  holes  (the  so-called  “multiple-site-damage” 
(MSD)),  are  examined  in  thorough  detail.  The  effects  of  various  non-dimensional  design  parameters  on  the 
reduction  in  the  stress-intensity  factors  near  the  crack-tip  are  determined,  and  are  presented  in  the  form  of 
design  charts.  Both  analytical  and  numerical  methods  are  employed  in  this  study.  In  the  analytical  method, 
the  cracked  metallic  plate  was  considered  to  he  infinitely  large  and  the  composite  patch  was  modeled  as  a 
long  orthotropic  strip  of  finite  height  (in  the  direction  perpendicular  to  the  crack  axis).  Next,  by  using  the 
FEAM,  a  more  general  analysis  capability  that  can  treat  arbitrary  shapes  of  the  cracked  metallic  sheet  as 
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Figure  3.1 :  Schematic  of  a  cracked  metallic  sheet  repaired  with  a  composite  patch 


well  as  the  composite  patches,  is  developed.  This  general  FEAM  is  applied  to:  (i)  composite  patch  repairs 
of  cracks  emanating  from  loaded  fastener  holes  (the  MSD  problem);  (ii)  composite  patch  repairs  of  semi¬ 
elliptical  shaped  surface  flaws  in  thick-section  plates;  and  (iii)  composite  patch  repairs  of  quarter-elliptical 
shaped  surface  flaws  emanating  from  fastener  holes.  The  problem  (i)  is  two-dimensional  in  nature,  while 
problems  (ii)  and  (iii)  are  fully  three-dimensional.  In  all  these  cases,  the  effects  of  various  design  parameters 
on  the  reduction  of  the  crack-tip  (front)  stress-intensity,  due  to  repair,  are  fully  discussed. 

Repair  of  center  or  edge  cracks  in  unstijfened  panels 

Due  to  the  fact  that,  in  general,  an  analytical  (as  opposed  to  purely  ni;  erical)  approach  is  more  conveniently 
employed  to  examine  the  effects  of  each  design  parameter  on  the  mechanical  behavior  of  a  repaired  crack, 
a  model  problem  of  a  repaired  crack  is  first  analyzed.  Fig.  3.2  depicts  the  problem  of  an  infinite  isotropic 
(metal)  sheet  containing  a  crack  of  length  2a,  repaired  by  an  adhesively  bonded  infinitely  long  orthotropic 
(composite)  patch  of  height  2H  (height  is  measured  normal  to  the  crack-axis).  It  is  assumed  that  the  crack 
lies  along  the  it-axis  and  that  the  principal  material  directions  of  the  patch  coincide  with  the  geometric  it  and 
y  axes.  Let  the  thickness  of  the  metal  plate,  the  composite  patch  and  the  adhesive  layer,  be,  respectively,  hs, 
hp  and  ha- 

The  metal  sheet  is  subjected  to  a  hoop  stress  of  Oo  tis  in  Fig.  3.1.  This  mode  I  symmetric  problem  can 
be  solved  by  a  superposition  of  the  two  following  problems: 

Problem  A:  Determine  the  stress  field  in  an  isotropic  metal  sheet  without  a  crack,  to  which  is  bonded  an 
orthotropic  patch  as  in  Fig.  3.2a. 

Problem  B:  When  the  tractions  at  the  location  of  the  crack,  as  determined  from  the  solution  of  Problem 
A  above,  are  erased,  determine  the  stress-intensity  factors  at  the  crack-tips.  This  problem  is  sketched  in 
Fig.  3.2b. 

The  relevant  material  and  geometric  parameters  from  Problems  A  and  B  are  as  follows: 

1.  The  Young's  modulus  Es,  and  Poisson's  ratio  v,.  for  the  isotropic  sheet. 


a 


b 


Figure  3.2:  Linear  superposition  of  problem  A  and  problem  B  to  solve  the  problem  in  Fig.  3.1 


2.  Thickness  of  the  isotropic  sheet. 

3.  The  Young's  modulus  and  Ey  of  the  orthotropic  patch;  the  Poisson  ratio  Vy  and  the  shear  modulus 
Gxy  of  the  orthotropic  patch. 

4.  The  height  H ;  and  thickness  hp  of  the  orthotropic  patch;  and 

5.  The  thickness  ha,  and  the  shear  modulus  Ga  of  the  adhesive  layer. 

By  using  the  Buckingham' s  0  theorem,  the  grouping  of  non-dimensional  parameters  that  arise  from  the 
above  geometric  and  material  data  is  determined  to  be: 


It  is  the  objective  of  the  present  analysis  to  determine  the  stress-intensity  factor  as  a  function  of  the 
above  8  non-dimensional  parameters.  This  ensemble  of  curves  would  then  lead  to  a  design  procedure  to 
select  the  appropriate  composite  patch,  and  adhesive,  for  a  given  crack  in  a  metallic  sheet.  This  is  the  first 
comprehensive  treatment  of  this  problem  in  tbe  literature. 

Problem  A  (Fig.  3.2a) 

Since  the  base  plate  and  the  patch  have  infinite  lengths  in  the  jc-direction  and  since  there  is  no  crack,  we 
can  treat  this  problem  as  a  one-dimensional  patch  problem.  Here  it  is  assumed  that  the  shear  stresses  in  the 
adhesive  layer  are  only  in  the  y-direction  (i.e.,t^2  =  0),  the  base  plate  and  the  patch  are  in  states  of  plane 
stress,  and  that  the  shear  stresses  in  the  adhesive  layer  are  applied  on  the  base  plate  and  the  patch  as  body 
forces  (Fig.  3.3). 

The  equilibrium  equations  for  a  metal  sheet  and  the  composite  patch  are,  respectively: 
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Figure  3.3;  Basic  problems  needed  in  solving  the  problems  of  Fig.  3.2b 


dOyy  Xy; 

dy  hp 


(3.3) 


Here  and  are  the  normal  stresses  in  the  y  direction  in  the  base  sheet  and  in  the  patch  respectively, 
and  is  the  shear  stress  in  the  adhesive  layer. 

The  stresS'Strain  relations  of  the  base  sheet  and  the  patch  are  given  by: 
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Here  Es  and  are  Young’s  modulus  and  Poisson’ s  ratio  of  a  base  shcei  respectively;  and  Ey  and  Vy  are 
the  moduli  and  Poisson’s  ratios  of  an  orthotropic  patch.  Among  these  variables,  the  relation  Vx/E>x  = 
is  satisfied. 

As  a  compatibility  condition  between  the  metal  sheet  and  the  composite  patch,  the  following  relations 
proposed  by  Mitchell,  Wooley,  and  Chivirut  (1975),  Jones  and  Callinan  (1979)  were  used 

Tv.  =  -j-  (3.8) 

Here 


Ga  8G\.  8C,, ' 


(3.9) 


Vj  and  vP  are  the  y  displacements  in  the  sheet  and  the  plate  respectively,  and  G^,  Gs  and  Gp  are  the  shear 
moduli  of  the  adhesive,  the  base  sheet  and  the  patch,  respectively.  We  assume  that  Gp  has  the  same  value 
as  Gxy  in  the  patch  plate.  From  the  elaborate  numerical  results  presented  later,  it  will  be  seen  that  this 
assumption  has  little  effect  on  the  stress  intensity  factors. 

By  using  Eq.  3.8  and  the  following  boundary  conditions: 
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(fyy  =  Go  at  V  =  //,  G^y  =  0  at  y  =  H, 
we  can  solve  Eqs.  3.2  and  3.3  easily.  The  stresses  are  obtained  as 


(3.10) 


^Jl=^^^(co,MH-ooshAy), 


(i-v?)<7o 
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(3.12) 

(3.13) 


Here,  to  simplify  the  problem,  the  additional  assumptions  =  0  and  =  Owere  used.  A^  is  defined  as 


^  =F 


1  /1-vf 


Echo 
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F  h 


(3.14) 


From  Eqs.  (3.11-3.13),  we  can  notice  that  and  have  their  maximum  values  aty  =  H,  and  decrease 
exponentially  as  {H  —y)  becomes  larger  than  zero.  On  the  other  hand,  G^  is  zero  aty  =  H,  but  it  increases 
exponentially  as  {H  —  y)  becomes  larger  than  zero. 

When  cosh/1//  1 ,  the  stresses  near  y  =  0  can  be  expressed  as 


G'yy  =  aO 


hz 


(1-v?) 

FEshsA^ 
FE,hpA^^^' 


cto, 


=  0. 


(3.15) 

(3.16) 

(3.17) 


These  stresses  are  the  same  as  those  when  the  patch  is  considered  to  be  of  an  infinite  height.  Thus,  when 
the  value  of  cosh  AH  is  much  greater  than  1 ,  we  can  assume  that  the  patch  has  an  infinite  height,  because  this 
assumption  gives  nearly  the  same  values  for  the  stress  intensity  factors  at  the  cracktip,  and  induced  stresses 
near  y  =  0. 


Problem  B  (Fig.  3.2b) 

In  order  to  solve  this  problem,  a  superposition  of  the  basic  problems  shown  in  Fig.  3.3a,  b  and  c  may  be 
used. 

Fig.  3.3a  is  the  case  where  uniform  surface  pressure  po  is  applied  on  the  crack  surface  in  a  metal  sheet 
of  infinite  size.  Its  solution  is  well  known  [Gladwell  and  England  (1977)].  Let  the  displacement  field  of 
this  problem  be  =  f\  (x,y)  ,\F^  =  f2  (x,y).  Here  is  the  displacement  in  the  x  direction  and  v*  is  the 
displacement  in  the  y  direction,  in  the  sheet. 

Fig.  3.3b  is  the  problem  wherein  point  loads  X,Y  are  applied  symmetrically  at  the  points  (A:o,yo), 
(-xo.yo),  (xo, -yo)  and  (-xq,— yo)  on  an  infinite  isotropic  plate  with  a  central  crack.  Let  this  solution 
be  and  \F^  (x.  y).  The  solutions  can  be  found  in  Erdogan  and  Arin  (1972)  and  Roderick  (1980).  By  using 
the  solution  in  Roderick  (1980),  the  displacements  at  a  point  (x,y)  for  this  problem  can  be  expressed  as: 
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Figure  3.4:  Green's  function  solution  for  a  finite  orthotropic  plate,  through  superposition 


=  Hi\{x,y\xo,yo)X +  H\2{x.y,xo,yo)Y  (3.18) 

=  H2\{x,y;xo,yo)X  +  H22{x,y,xo,yo)y  (3.19) 

The  problem  of  Fig.  3.3c  is  where  the  point  loads  X,  Y  are  applied  symmetrically  at  the  points  (xo,yo). 
(— a:o,yo)»  (^Oi  — yo)  and  {—xq,  — yo)  on  an  orthotropic  patch  of  height  2H.  The  problem  of  Fig.  3.3c  can 
again  be  solved  through  a  linear  superposition  of  two  other  problems  as  in  Fig.  3.4b  and  c.  Fig.  3.4b  is  of  the 
problem  of  an  infinite  orthotropic  sheet;  and  1 6c  is  of  a  finite  height  composite  patch  with  residual  tractions, 
from  the  solution  of  the  problem  in  Fig.  3.4b,  acting  on  the  edges  y  =  ±/f.  The  displacement  field  for  the 
problem  in  Fig.  3.4b  can  be  found  in  Leknitskiic  (1968).  Now  we  consider  the  problem  of  Fig.  3.4c.  Here, 

-p{x)  and  -q{x)  are  the  normal  and  shear  stresses  acting  on  the  y  —  H  plane.  Since,  for  the  problem  in 

Fig.  3.4c,  the  displacements  uP^  (x,y)  and  vP‘^  (x,y)  are  odd  and  even  functions,  respectively,  with  respect  to 
X,  we  may  define: 


uP^  ix.y)  =  -  I  <j)(a.  vlsinavz/a.  (3.20) 

n  Jo 

vP‘-' {x.v)  =  -  [  V)  cos(.',  ,,  a,  (3.21) 

7t  Jo 

where  (|)  (a,y)  and  X|/(a,y)  are  the  Fourier  transforms  of  iiP‘'  [x.y]  and  v^‘^(x,y)  respectively.  Since  there  are 
no  further  body  forces,  the  equilibrium  equations  for  the  patch  in  Fig.  3.4c  are  written  in  terms  of  and 
vP'^  for  this  component  problem  of  Fig.  3.4c,  as: 


\l-VxVyJ  axdy  \l-VvVv/  dv-  \  dxdy  dx^  J 

and 


f  E,  \d\P\  (  Vv£v  _ 

Vl-V,vj  3x2  a.vdv  ''"A  3y2  ^dxdy)~^' 

Substituting  Eq.3.20-3.21  into  Eq.3.22-3.23,  we  obtain: 


(3.22) 


(3.23) 


0, 


(3.24) 
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The  general  solutions  of  Eqs.  3.24-3.25  are: 


(j)  =  ei/4i  coshPiy-l-g2^2Cosh(32y-l-ei/i3sinhPi}’-fe2>i4SinhP2)’,  (3.26) 

r|/  =  A|  sinhpi>'-|-A2  sinhP2y  +  ^3CoshPi>'-fy44cbshP2y7  (3.27) 


where 


_  OtP*  [v,r£^v  +  Catv  ( 1  '^ArVy)]  i  _  ^ 

G„pJ(l-v,v,)-a2E, 

and  Pi  and  P2  are  roots  of 


(3.28) 


P^  -  a"  (^  -  2v;,^  p2  +  =  0  (3.29) 

Considering  the  symmetry  conditions,  {x,y)  =  (x,  -y),  and  (x,y)  =  (x,  -y),  we  see  that 

the  constants /43  and  A4  must  be  zero  in  Eqs.  3.26-3.27.  Using  Eqs.  3.26-3.27  in  Eq.  3.20-3.21,  we  obtain 
the  displacements  uP‘^  and  From  these,  the  corresponding  stresses  can  be  determined  as: 


where 


2  f°° 

=  -J  [(aMi^i -f  M2Pi)Ai  coshPiy-H  (aA/ie2  +  A^2P2)>i2Coshp2y]cosax6fa 
2  /■“ 

o^y  =  -  [(aM2ei -1-A/3Pi)Ai coshPiy-H (aA/2e2  +  A^3p2)'^2CoshP2y]cosctxda 

7C  J  0 
2 

^xy  =  -  [G^.(eipi -a)>lisinhpiy  +  G^y(e2P2-a)/42sinhP2y]sin(Xxf/a, 

J  0 


A/,  =  — M,  =  M3  =  — ^ 

1  -  VjcVy  ■  1  -  GVyVv  1  -  VrVv  1  -  VyVv 


If  P  [a)  and  g  (a)  are  the  Fourier  transforms  of  the  tractions  p  (x)  and  q  (x)  in  Fig.  3.4c,  such  that 


(3.30) 


p{x)  = 

2 

K  . 

2 

/*oo 

1  P{a)cosaxdx, 

Jo 

roo 

(3.31) 

q{x)  = 

K  . 

1  Q  (a)  sin  oucdx, 

Jo 

(3.32) 

then  the  boundary  conditions  at  ayy  (jc,  H)  =  +P  [x]  and  (x,  H)  =  +q  (x)  result  in  the  following  equations: 

{aM2e\  +A/3pi)  (coshPi//)Ai  +  (aAf2e2  +  ^f3P2)  (coshp2/7)A2  =  P(a)  (3.33) 


Gxv(^iPi  -a)  (sinhp|//)y4i -|-G;cy  (g2P2 -oc)  (sinhP27/)/42  =  2(a)  (3.34) 

from  which  /li  and  A2  are  solved  for.  Thus,  the  displacement  solution  (x,y)  and  (x,y)  is  obtained  for 
the  problem  in  Fig.  3.4c. 
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Let  the  displacement  solution  for  the  patch  problem  of  Fig.  3.3c,  obtained  as  a  superposition  of  problems 
in  Fig.  3.4b  and  c,  be  written  as: 

=  ~  j  [A'ii(x.y;xo,yo)'t;cz(xo,yo)  +  A'i2(x,y;xo,yo)%(-^o,>'o)]^^^of5?yo  (3.35) 

J  D 

and 

=  ]^  [  [K2i{x,r,xo,yo)Txz{xo,yo)+K22{x,y,xo,yo)Xyz{xo,yo)]dxodyQ  (3.36) 
J  D 

wherein,  the  adhesive  shear  stresses  (x^z/hp)  and  (Xyz/hp)  are  treated  as  body  forces  per  unit  volume,  acting 
on  the  composite  patch,  as  in  Fig.  3.3c. 

The  sheet  displacements  from  problems  3.3A  and  3.3B  may  be  written  as: 


{x,y)  +  M''®  (x.y) 

/i  [  [^11  (j^.y;xo,yo)  Xxz  (j:o,.vo)  +  (x,>’;xo,yo)  V  ixo,yQ)]dxodyQ  (3.37) 


and 


(-^ly)  =  /2  (^,y)  ~]^  [^21  (■*,y;-ro,.yo)  Xxz  (xo,yo)  +  H22  (j:,y;yo,yo)  %  (-yo.yo)]  dx^dyo  (3.38) 

where  D,  the  domain  of  integration,  is  the  size  of  the  composite  patch.  The  body  forces  on  the  sheet  are 
{—Xxz/hs)  and  {—Xyz/hs)  in  the  x  and  y  directions,  respectively. 

The  compatibility  condition  between  the  sheet  and  the  patch  is  expressed  as 


and  Xv-  = 


(3.39) 


By  substituting Eqs.  3.35-3.36  and  Eqs.  3.37-3.38  into  Eq.  3.39.  one  v  oiains  the  following  integral  equations 
for  Xxz  %• 


PXxz  (jc^y)  +  [^1 1  (x,y;xo,yo)  x^,-,  (xo.yo)  ^'12  (.v.y:.V(i.y())  x,.,  (xo,yo)]^/xoc?yo  =  f\  (x,y)  (3.40) 

and 

Phz  {x,y)  +  (yfyi^o^yo)  Tv-  (xo-.vo)  -b  a-::  ':.v.y:.v,).yo) Xy, (xo,yo)]^^xo^/yo  =  fi (-y,y) 


where 


^11 


('r+'r) 

\  hg  hp  J 


ki2  = 


V  h,  ^  hp  J 


/C2t  = 


[^s  hpj- 


(3.41) 


In  integral  Eqs.  3.40  and  3.1.1,  the  domain  of  integration  D  is  the  size  of  the  orthotropic  patch,  which 
is  finite  in  the  y  direction,  but  «>  in  the  x  direction.  For  numerical  purposes,  however,  the  domain  D  in 
the  X  direction  is  also  truncated  to  some  finite  value.  In  order  to  solve  Eqs.  3.40  and  3.1.1,  the  domain 


W/a  =  4 


Crack 


Figure  3.5:  A  typical  finite  element  mesh  of  the  adhesive  layer 


D  is  discretized  into  a  number  of  small  sub-elements.  In  each  subelement,  and  Xyz  are  assumed  to  be 
constant.  Thus  a  system  of  linear  equations,  for  discrete  values  of  and  Xy^,  is  obtained.  By  solving  these, 
the  shear  stresses  Xxz  and  are  determined.  Since  the  stress-state  on  the  sheet,  as  in  Fig.  3.3a  and  b,  is  thus 
completely  known,  the  /c-factors  can  be  evaluated  from  the  analytical  solutions  for  problems  3.3A  and  B. 

One  example  of  descritization  of  D  is  shown  in  Fig.  3.5,  for  the  case  when  (H/a)  =  4.  By  varying  the 
size  of  this  D,  it  was  determined  that  the  obtained  A:-factors  are  insensitive  to  the  size  of  D,  as  long  as  the 
dimensions  of  D  are  about  five  times  as  large  as  a.  Note  also  the  finer  mesh  near  the  crack,  in  Fig.  3.5,  as 
the  adhesive  shear  stresses  are  expected  to  be  maximum  near  the  crack  itself. 

Behavior  of  patched  cracks  in  unstiffened  panels 
Attention  is  focused  in  this  study,  on  the  variation  of: 

1 .  non-dimensional  stress  intensity,  i.e.,  Ki/ooy^ 

2.  non-dimensional  adhesive  shear,  ^  and  — 

<^0  0() 

3.  non-dimensional  patch  stress,  [Oyyhp]  /  (00^^) 

as  a  function  of  the  non-dimensional  parameters: 


From  the  computed  results,  it  was  found  that  the  more  significant  parameters  are: 

1 .  ■  which  characterizes  the  adhesive  flexibility.  The  “smaller”  this  number  is,  the  “stiffer”  is  the 

adhesive. 
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■  which  characterizes  the  patch  stiffness  in  the  direction  normal  to  the  crack  axis. 
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Figure  3.6:  (a)  Reduction  in  normalized  /:-factor  due  to  patching,  as  a  function  of  the  crack  length;  (b) 
reduction  in  the  actual  A:-factor  (for  =  1)  due  to  patching  as  a  function  of  the  crack  length 


The  “larger”  this  number,  the  “stiffer”  is  the  patch  in  the  direction  normal  to  the  crack-axis. 

Fig.  3.6a  shows  the  variation  of  the  normalized  stress-intensity  factor  as  the  crack-length  increases,  for 
the  case  of  a  very  large  composite  patch,  {H/a)  — >  {hpEy/hsGs)  —  2;  and  {GsF/hs)  =  3.0.  Note  that, 
in  the  absence  of  a  composite  patch,  the  normalized  A:-factor,  i.e.,  \Kj/a\/%S[  =  1.0.  Thus,  the  vertical 
distance  between  the  horizontal  straight  line  [normalized  /c-factor  =  1.0]  for  the  unpatched  crack,  and  the 
curve  for  the  repaired  crack,  indicates  the  reduction  in  the  normalized  k-factor  due  to  patching.  Fig.  3.6a 
clearly  shows  that  more  and  more  reduction  in  the  normalized  /:-factor  results  as  the  crack-length  increases. 

Also  it  is  seen  from  Fig.  3.6a  that,  as  [a/hs)  <^\,  the  normalized  ^-factor  converges  to  that  in  a  center- 
cracked  sheet  with  the  crack-face  pressure  of  pQ  (see  Fig.  3.3a),  which  is  smaller  than  the  applied  far-field 
stress  ao  in  the  sheet  (see  Fig.  3.2a).  Thus,  when  the  crack  length  in  the  metal  sheet  is  very  small,  the 
reduction  in  k  -factor  is  achieved  by  the  reduction  in  the  crack-face  pressure,  po,  in  the  metal  sheet,  due  to 
the  pressure  of  the  composite  patch.  From  Eq.  3.1 1  it  is  seen  that  ;n  Fig.  3.2b  is  given  by 

PO  =  o;,  fr  =  0)  =  ao  -  (coshAH  - 1)  (3.42) 

Substituting  for  from  Eq.  3.14  into  3.42,  it  is  seen  that  pQ  is  independent  of  the  adhesive  flexibility 
parameter  “FV  for  large  values  of  H  when  cosh  AH  >  1 .  Thus,  for  short  cracks,  the  reduction  in  the 
normalized  A:-factor  is  unaltered  by  the  adhesive  flexibility.  As  the  crack  becomes  longer,  the  adhesive  shear 
stresses  acting  on  the  metal  sheet  have  more  dominant  effect  on  the  crack-tip  A:-factor,  and  the  reduction  in 
the  A:-factor  depends  on  the  adhesive  flexibility  parameter  F. 

Fig.  3.6b  shows  the  variation  of  the  actual  ^-factor  as  the  crack  grows,  with  and  without  the  presence 
of  a  composite  patch.  It  is  seen  that  without  the  patch,  the  ^-factor  increases  as  y/a  as  expected,  but  in 
the  presence  of  a  patch,  the  /:-factor  levels  off  to  a  constant  value  as  the  crack  grows.  This  result  has  also 
been  noted  earlier  by  Rose  (1981)  and  Atluri  and  Kathiresan  (1978).  It  should  be  noted,  however,  that  in 
the  present  analysis,  the  adhesive  bond  between  the  metal  sheet  and  the  composite  patch  is  assumed  to  be 
perfect.  However,  the  shear  stress  in  the  adhesive  is  maximum  near  the  center  of  the  crack,  where  the  crack¬ 
opening  displacement  v*  in  the  sheet  is  maximum.  Thus,  for  all  adhesives  whose  maximum  shear-strength 
is  a  finite  value,  a  disbond  can  be  expected  to  occur  between  the  metal  sheet  and  the  composite  patch  near 
the  center  of  the  crack  at  y  =  0,  as  well  as  possibly  at  y  —  H.  When  this  disbonding  is  modeled  properly,  its 
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Patch  stiffness  parameter  hpEy/hsGs 


Figure  3.7:  Reduction  in  normalized  ^-factor  for  a  given  crack  length,  as  a  function  of  the  patch  stiffness 
and  adhesive  flexibility 


effect  is  to  increase  the  actual  /c-factor  of  the  repaired  crack,  as  the  crack  length  increases.  This  situation  of 
the  damage  tolerance  of  the  bond  is  being  currently  quantified. 

Fig.  3.7  shows  the  effects  of  the  “patch-stiffness”  parameter  and  the  “adhesive-flexibility”  parameter  on 
the  normalized  stress-intensity  factor  of  the  patched  crack,  for  a  given  crack  length  (a/hs)  =  5.0,  and  a  given 
patch-height,  {H /a)  =  4.0. 

For  a  fixed  “adhesive-flexibility”  value,  the  stress-intensity  factor  (S.I.F.)  of  the  patched  crack  decreases 
rapidly  from  its  “unpatched”  value,  as  the  patch-stiffness  increases.  This  S.I.F.  reduction  is,  however,  smaller 
as  the  adhesive  flexibility  increases. 

From  Fig.  3.7,  it  is  evident  that  the  S.I.F.  reduction  due  to  a  patch  is  maximum  when  one  uses  a  very  stiff 
composite  patch,  as  well  as  a  very  stiff  adhesive  layer. 

However,  a  stiff  adhesive  layer  has  a  deleterious  effect  on  the  integrity  of  the  bond  itself  as  the  adhesive 
shear  stresses  increase  in  such  a  case,  as  seen  from  Fig.  3.8  and  3.9. 

Fig.  3.8  shows  the  non-dimensional  shear  stress  {XyJ Oq)  in  the  adhesive  layer,  at  the  center  of  the  crack, 
i.e.,  x  =  0  and  y  =  0;  while  Fig.  3.9  shows  the  corresponding  value  at  the  end  of  the  patch,  i.e.,  x  =  0  and 
y  —  H.  It  is  noted  that  in  all  cases,  Xxz  C  Tyj  and  hence  Xxz  ‘S  not  shown.  It  is  seen  that  Xy^  at  y  =  0  arises 
only  from  the  problem  in  Fig.  3.2b.  Since  ty,  =  (v'  -  v^)  /F  and  since  (v'  —  v^)  is  maximum  at  x  =  0  and 
y  =  0  in  Fig.  3.2b,  the  maximum  value  for  the  problem  in  Fig.  3.2b  also  occurs  at  x  =  0  and  y  =  0.  Also, 
the  stress  Xyz  at  x  =  0  and  y  =  //  for  the  problem  in  Fig.  3.1  are  nearly  the  same  as  that  for  the  problem  in 
Fig.  3.2a,  and  given  by  Eq.  3.13  when  y  =  H. 

From  Fig.  3.8  it  is  seen  that  Xyz  at  x  =  0  and  y  —  Q  increases  rapidly  as  the  adhesive  becomes  stiffer  (i.e., 
as  F  decreases),  especially  for  moderate  values  of  patch  stiffness.  For  fixed  values  of  F,  (Xyz/oo)  at  x  =  0 
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Patch  stiffness  parameter 


Figure  3.8:  Variation  of  at  j:  =  0,  j  =  0  as  a  function  of  the  patch-stiffness  and  adhesive  flexibility 


and  y  =  0  decreases  as  the  patch  stiffness  increases.  For  large  values  of  adhesive  flexibility,  the  maximum 
shear  stress  is  nearly  constant  regardless  of  the  value  of  the  patch  stiffness. 

From  Fig.  3.9  it  is  seen  that  x^^at  x  =  0  and  y  =  H  also  increase  is  the  adhesive  layer  becomes  stiffen 
However,  for  fixed  values  of  F,  this  shear  stress  increases  as  the  pat^n  increases. 

Thus,  we  note  that  when  a  patch  with  a  low  stiffness  {  [hpEp]  /  [hsGs)  <  3.0)  is  used,  the  value  of  Xy^  at 
X  =  y  =  0  is  larger  than  the  value  of  x^^  at  a:  =  0  and  y  =  H .  As  the  patch  becomes  stiffer,  x^^  at  a:  =  0  =  y 
decreases,  while  Xyz  at  =  0  and  y  =  H  increases.  Thus  as  the  patch  becomes  very  stiff,  Xy^{x  =  0,y  —  H)  > 
yz{x  =  0,y=0). 

Thus,  in  general,  depending  of  the  values  of  the  patch  stiffness,  the  likelihood  of  a  disbond  between  the 
metal  plate  and  the  patch  exists  both  at  the  locations  (.v  =  0.  y  =  0)  and  {x  =  0,y  =  H). 

Fig.  3.10  shows  the  effect  of  the  patch-stiffness  and  adhesive  flexibility  on  the  values  of  {Cyyhp)  / 

The  value  of  K/  ao)  changes,  on  the  other  hand,  according  to  the  value  of  {hs/hp).  For  fixed  value  of 
adhesive  flexibility,  [O^yyhp  )  /  (cohs)  increases  as  the  patch-stiffness  increases.  For  fixed  patch-stiffness, 
[a^h  p)  /  i^ohs)  increases  as  the  adhesive  stiffness  increases. 

So  far,  we  have  assumed  that  the  metal  sheet  is  infinite,  and  the  composite  patch  is  of  a  finite  height  but 
infinitely  long.  In  reality  however,  the  base  sheet  as  well  as  the  patch  are  finite  in  size;  moreover,  the  patch 
can  be  arbitrary  in  shape.  To  account  for  these  size  and  shape  effects,  one  may  use  the  FEAM  to  determine 
the  displacement  fields  in  each  of  the  component  problems  such  as  in  Fig.  3.3a,  b,  c  wherein  each  problem 
may  have  a  finite,  arbitrary,  shape.  The  integral  equations  for  tve  in  these  finite  arbitrary  shapes  may 

be  formulated  in  much  the  same  way  as  described  earlier. 
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Patch  stiffness  parameter  h^Ey^/hsC^ 

Figure  3.9:  Variation  of  Xy^  si  x  =  0,  y  =  H  as  a  function  of  the  patch-stiffness  and  adhesive  flexibility 


3,1.2  Repair  of  cracks  near  loaded  fastener  holes  in  a  lap  joint  (multiple-site  damage) 

Consider  the  model  problem  of  repairing  of  symmetrical  cracks  emanating  from  a  loaded  fastener  hole,  as 
shown  in  Fig.  3.1 1. 

A  circular  hole  of  radius  R  is  located  at  the  center  of  a  rectangular  plate  of  width  IWs,  height  2//^  and 
thickness  hs.  Two  radial  cracks  of  equal  lengths  [a  -  R)  each  emanate  from  the  hole,  along  the  it  axis.  To 
repair  these  cracks,  a  rectangular  orthotropic  patch  of  width  2Wp,  height  2Hp  and  thickness  hp  is  assumed 
to  be  bonded  onto  the  metal  sheet.  The  thickness  of  the  adhesive  is  ha.  The  metal  sheet  is  assumed  to  be 
subject  to  a  far-held  hoop  stress  Gq  as  well  as  the  pin-loading  p  (0)  of  the  type: 

P(0)  =  ^^^|sin0|  (3.43) 

as  Fig.  3.12a.  Since  the  applied  loading  (ctoand  p(0))  is  not  symmetric  in  the  y  direction,  the  crack  may 
exhibit  mode  II  behavior  also.  However,  this  mode  II  k-factor  may  be  expected  to  be  much  smaller  than  the 
mode  I  component,  and  will  be  ignored  hence  forth. 

The  problem  of  Fig.  3. 1 2a  can  be  solved  by  a  superposition  of  the  two  symmetric  problems  of  Fig.  3. 1 2b 
and  c.  If  the  crack-tip  /:-factor  in  Fig.  3.12b  is  Kjb  and  that  in  Fig.  3.12c  is  Kic,  respectively,  then  the  k 
-factor  of  Fig.  3. 12a  may  be  written  as; 


=  2  +  ^ic)  ■ 

Note  that  both  the  problems  of  Fig.  3. 1 2b  and  c  involve  finite  geometries. 


(3.44) 
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Patch  stiffness  parameter  hpEy/h^Gs 


Figure  3.10:  Variation  of  at  X  =  0,  >>  =  A/  as  a  function  of  the  patch-stiffness  and  adhesive  flexibility 


The  problems  in  Fig.  3.12b  and  c,  can,  in  turn,  be  solved  by  a  linear  superposition  of  the  basic  problems 
in  Fig.  3.13:  (A)  the  problem  of  a  finite  sized  isotropic  metal  sheet  with  radial  cracks  emanating  from  an 
unloaded  hole,  the  sheet  being  subject  to  far-field  hoop  stress  Oo;  ( B  ihe  problem  of  a  finite-sized  isotropic 
metal  sheet  with  radial  cracks  emanating  from  a  hole,  the  hole  surtae .  being  subject  to  symmetric  load  p  (0) 
(Fig.  3.13b);  (C)  the  problem  of  a  finite-sized  isotropic  metal  sheet  with  cracks  emanating  from  a  hole,  the 
plate  being  subjected  to  shear  tractions  Txz  and  Xy,  (from  the  adhesive  layer)  on  its  surface  (Fig.  3.13c);  and 
(D)  the  problem  of  a  finite-sized  orthotropic  composite  patch,  subjected  to  shear  tractions  ixz  and  lyz  (from 
the  adhesive  layer)  on  its  surface  (Fig.  3.13d). 

In  order  to  solve  the  above  four  basic  problems,  the  FE.AM  is  employed.  In  the  case  of  the  problem  in 
Fig.  3.13a,  the  FEAM  procedure  is  as  follows: 

1.  First  consider  the  analytical  solution  for  the  problem  of  a  central  crack  in  an  infinite  sheet,  the  crack- 
faces  being  subjected  to  arbitrary  tractions.  In  the  present  problem  of  cracks  emanating  from  a  hole,  as 
in  Fig.  3.13a,  in  the  FEAM  procedure,  residual  tractions  would  arise  on  the  crack-faces,  |a:|  >  {R  —  a) 
(for  a  traction-free  hole),  and  further,  there  will  be  steep  stress  gradient  near  |jc|  =  R  because  of  the 
stress-concentration  near  the  hole.  Thus,  in  order  to  use  the  Gladwell  and  England  (1977)  solution, 
tractions  on  the  fictitious  crack-face  |x|  <  a  will  have  to  be  interpolated  with  Chebyschev  polynomials 
which  is  not  desirable,  since  the  residual  tractions  on  the  fictitious  crack-faces  |x|  <  are  zero,  while 
those  on  |x:|  >  {Ra)  have  a  steep  gradient.  Due  to  this  fact,  an  alternate  analytical  solution  for  a  central 
crack  in  an  infinite  sheet,  the  crack  faces  being  subject  to  equal  and  opposite  point-loads  at  an  arbitrary 
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Figure  3.13:  Basic  problems  to  be  solved  in  solving  the  problems  of  Fig.  3.12 


location  x,  is  considered. 

Without  loss  of  generality,  consider  the  problem  when  a  crack  of  length  2a  exists  in  an  infinite 
isotropic  plate,  and  two  point  forces  of  magnitude  Y  are  applied  at  points  x  —  id  on  the  upper  crack 
surface,  and  two  point  forces  —  F  are  applied  at  x  =  ±d  on  the  lower  surface.  Following  Muskhel- 
ishvili  (1953),  the  complex  stress-functions  for  this  problem  are: 


0(z)  =  Q(z) 


/  1  1  \ 
IniVz^  -  a~  \s:-l-d  ^  z-d) 


(3.45) 


<pU)  =  <«(z)  =  ^ 


Here 


2-\/d~  —  a~Vz~  —  a~  +  2dz  -  2cr  2\/ d^  ~  cd-y/z}-  —  —  2dz  —  2a~ 

log - ; - 1- log 


z-d 


z+d 


(3.46) 


0(z)  =  (p'(z),  Q(z)  =  (o'(z). 

The  stress  intensity  factor  is  given  by: 


K,  =  ^( 

\/ tCfl  \  \/  ^/  —  cl  a  -\-  d  / 

The  stresses  and  displacements  are  obtained  from  the  following  relations: 


(3.47) 


1 


0(:,i-0(:) 


<yyy  -  i\xy  =  O  ( C  j  -r  )  O'  (z)  4-  Q  (z) 


(3.48) 

(3.49) 


and 


2G(i/  +  /v)  =  K(p(z)  -co(z)  -  (z-z)O(z) 


(3.50) 
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By  using  the  above  solutions  as  Green's  functions,  one  can  obtain  the  stress-intensity  factors,  and  the 
stress  and  displacement  fields,  when  arbitrary  tractions  are  applied  on  the  crack  surface. 

The  main  steps  of  the  FEAM  for  solving  the  problem  of  Fig.  3.13a  are: 

2.  Using  the  analytical  solution  for  the  problem  of  an  infinite  sheet  with  a  hole  (but  no  cracks)  subjected 
to  far-field  stress  Gq  (see  Muskhelishvili  (1953)),  obtain  the  residual  tractions  at  the  outer-boundaries 
of  the  finite  sheet  as  well  as  at  the  locations  of  the  crack,  |x|  >  {a  —  R). 

3.  To  create  a  traction-free  crack,  erase  the  stresses,  as  found  in  step  (2),  on  the  crack-faces  |x|  >  {a  —  R), 
using  the  analytical  solution  for  an  infinite  sheet  with  a  crack,  as  in  Eqs.  3.47-3.50.  Determine  the 
A:-factors,  and  the  displacement  solutions  using  the  Green’s  function  solution,  Eqs.  3.47  and  3.50. 

4.  Corresponding  to  the  solution  in  step  3,  determine  the  residual  tractions  at  the  surface  of  the  hole  as 
well  as  at  the  outer  boundaries  of  the  finite  plate. 

5.  In  order  to  satisfy  the  given  traction  b  e  at  the  outer  boundaries  of  the  finite  plate  as  well  as  at  the  hole 
surface,  reverse  the  residual  tractions  at  these  surfaces,  as  calculated  from  steps  (2)  and  (4).  Calculate 
the  equivalent  nodal  forces  on  these  surfaces,  at  the  chosen  finite  element  nodes. 

6.  Using  the  usual  finite  element  method,  obtain  the  stresses  at  the  location  of  the  crack,  and  the  dis¬ 
placement  field  everywhere,  corresponding  to  the  nodal  loads  in  step  (5). 

7.  Erase  the  crack-face  stresses  as  computed  from  step  (6),  by  repeating  step  (3). 

8.  Continue  the  iteration  until  the  residual  tractions  in  the  repeated  step  (6)  are  vanishingly  small. 

9.  By  summing  appropriate  contributions,  compute  the  total  stress-intensity  factor,  and  the  total  dis¬ 
placement  field  for  the  problem  in  Fig.  3.13a. 

The  procedures  for  the  solutions  of  the  problems  in  Fig.  3. 1 3b  and  c  are  entirely  analogous  to  that  in 
solving  the  problem  in  Fig.  3. 1 3a,  as  described  above,  except: 

10.  In  solving  the  problem  of  Fig.  3.13b,  the  analytical  solution  in  step  (2)  of  the  problem  of  Fig.  3.13a 
is  replaced  by  an  analytical  solution  for  an  infinite  sheet  with  a  hole  (but  no  cracks),  the  hole  surfaces 
being  subject  to  symmetric  pin-loading,  p  (0). 

1 1.  In  solving  the  problem  of  Fig.  3.13c,  the  analytical  solution  in  step  (2)  of  the  problem  of  Fig.  3.13a 
is  replaced  by  an  analytical  solution  for  an  infinite  sheet  with  a  hole  (but  no  cracks),  the  sheet  being 
subjected  to  body  forces  and  x^^  a  rising  from  the  adhesive  layer.  (This  analytical  solution  is 
obtained  using  the  complex  variable  information,  as  in  Muskhelishvili  (1953)). 

In  the  case  of  the  problem  of  Fig.  3.13d,  the  following  steps  are  used: 


18 


(a)  By  using  the  analytical  solution  [Leknitskiic  (1968)]  of  the  problem  of  an  infinite  orthotropic 
sheet  subjected  to  four  sets  of  point  loads  applied  symmetrically  at  (±a:o;  ±yo).  we  obtain  the 
residual  stresses  at  the  outer  boundaries  of  the  patch,  x  =  dz/Z^;  and  y  —  ±.Wp,  as  well  as  the 
displacement  fields. 

(b)  In  order  to  satisfy  the  traction-free  boundary  conditions  at  a:  =  ±Hp  and  y  =  ±.Wp,  reverse  the 
residual  tractions  computed  in  step  (i)  above,  and  compute  the  equivalent  nodal  forces  to  be 
applied  at  the  (coarse)  finite-element  mesh  of  the  composite  patch. 

(c)  Using  the  finite  element  method,  obtain  the  displacement  field  corresponding  to  the  loads  in  step 
(ii).  Add  the  displacement  fields  of  steps  (i)  and  (iii). 

Formulation  of  integral  equations  for  the  problem  of  repair  of  cracks  near  loaded  fastener  holes 

By  using  the  FEAM  outlined  in  Sec.  3.1.2,  we  obtain  the  displacement  fields  in  each  of  the  four  basic 
problems  sketched  in  Fig.  3.13a-d,  respectively.  Let  the  displacement  field  of  the  problem  in  Fig.  3.13a  be 
given  by  u^^  =  Fia  (x,y)  and  =  F2a  (•v,y);  and  that  for  problem  in  Fig.  3.13b  be  given  by  u^^  =  Fib  (jC)}’) 
and  =  F2b  Let  the  displacement  field  for  the  problem  in  Fig.  3.13c  be  given  by: 

=  -"T  /  [f'fi\ix,r,xo,yo)^xz{xo:yo)  +  hi\2  (•^,.v;A:o,yo)%(j:o,yo)]^?^o^^yo  (3.51) 

«5  JD 

and 


(jc,y)  =  -  ^  ^  [H2\  (jc,y;j:o,.vo)  Tv;  (xo,yo)  +]dxodyo 


(3.52) 


where  D  is  the  domain  of  the  finite  patch  (and  of  the  adhesive  layer).  Likewise,  the  displacement  field  in  the 
patch  is  given  by 

uP  =  u’’^  {x,y)  -  ^  f  [Kii(x,y,xo,yo)Xx:.{xo.yo)~  Kizix.,  ).yo)Tv::  (jco,yo)]<^JCo^yo  (3.53) 
hp  J  D 

and 


yP  =  yP^  (x,y)  =  7^  /  [Zf2i  (^,yuo..vo)Tv;  (-ro-.vo)  -  K22  ix.y;xo,yo)Tyz{xo,yo)]dxQdyo  (3.54) 
h-n  J D 


The  total  displacement  in  the  metal  sheet  (see  Fig.  3.1 3a-c )  is  given  by: 


1 


sC 


and 


(3.55) 


=  1  ( +  r'-'-  =  1  (f.,,  +  f2B)  +  (3.56) 

The  conditions  of  compatibility  between  the  metal  sheet  and  the  orthotropic  patch  are  given,  again,  by: 

id  — 


and  Tv-  = 


F 


(3.57) 
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Figure  3.14:  A  typical  finite  element  mesh  of  the  adhesive  layer 


where  F  is  the  adhesive  flexibility  as  defined  in  Eq.  3.9.  By  substituting  Eqs.  3.51  to  3.56  in  Eq.  3.57,  one 
obtains  the  integral  equations  for  and  iy^.  Once  these  integral  equations  are  solved  for  Xxz  and  the 
free  -body  diagram  of  the  stress-intensity  factors  near  the  crack-tips,  as  effected  by  the  repair  patch,  can  be 
computed. 

By  assuming  constant  values  of  Xxz  and  within  each  subdomain,  the  integral  equations  in  Eqs.  3.56 
and  3.57  are  solved  for.  One  example  of  the  mesh  used  for  solving  the  integral  equation  over  D,  for  the  case 
Wp/R  =  2  and  Hp/Wp  =  2.0  is  shown  in  Fig.  3.14. 

In  the  problems  of  Fig.  3.12,  the  ^-factor  of  the  problems  in  Fig.  3.1 2a-c  are  identified  as  Kia  ,  Kjb  and  Kjc 
respectively,  [Kja  =  ^  {Kib  +  F-jc))-  In  the  following,  the  results  for  the  normalized  ^-factors,  as  functions 
of  the  material  parameters  (GsF/hs),  {hpEy/hsGs)  and  the  geometrical  parameters  a/R,  R/hs,  W/R,  and 
H/W,  are  discussed  for  the  patched  and  unpatched  cases. 

Results  and  discussions 

In  order  to  examine  the  accuracy  of  the  presently  developed  FEAM,  first  the  problem  of  cracks  emanating 
from  a  loaded  fastener  hole,  without  the  presence  of  the  repair  patch,  is  solved  first.  The  present  results  are 
compared  with  those  of  Cartwright  and  Parker  (1972),  for  the  case  when  H /W  =  2,  and  W /R  =  4.0.  In  the 
present  FEAM,  65  eight-noded  isoparametric  elements  were  used,  and  the  CPU  time  required  to  solve  one 
crack  case  was  about  55  seconds  on  a  Micro  VAX  Station  11.  The  present  results  are  seen  from  Fig.  3.15  to 
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Figure  3.15:  Variation  of  normalized  ^-factor  as  a  function  of  the  crack-length  (unrepaired  case) 


agree  excellently  with  those  of  Cartwright  and  Parker  ( 1 972). 

Fig.  3.16a  shows  the  effect  of  patching  on  the  normalized  strc'-s-intensity  factor,  as  the  crack  length 
increases,  for  given  material  and  geometrical  parameters  of  the  patci  .  as  identified  in  the  inset  of  Fig.  3.16a. 
Comparing  Fig.  3.15  and  3.16a,  the  reduction  in  the  normalized  A-factor,  as  the  crack  grows,  due  to  the 
presence  of  the  composite  patch,  can  be  observed.  For  small  cracks  this  reduction  is  not  much,  but  as  [a/R) 
increases,  the  effects  of  patching  increases  as  well. 

Fig.  3.16b  shows  the  effect  of  patching  on  the  actual  stress-intensity  factor,  as  the  crack  length  increases, 
for  given  geometrical  and  material  parameters  as  identified  in  the  inset  of  Fig.  3.16b.  Once  again  it  is  seen 
that  even  for  a  repaired  crack  near  a  loaded  fastener  hole,  the  actual  k  -factor  becomes  a  constant  as  the 
crack-length  increases  under  the  patch.  This  appears  to  be  a  salient  feature  of  all  cases  when  the  crack  in  a 
fuselage  skin  is  repaired  with  a  stiff  composite  patch.  Fig,  3. 1 7  shows  the  effect  of  the  hole  radius  on  the 
normalized  k-factors  for  the  repaired  cracks. 

Fig.  3.18  shows  the  effect  of  the  adhesive  layer  flexibility  and  patch-stiffness  on  the  normalized  stress- 
intensity  near  the  crack-tip,  for  the  case  of  far-field  hoop-stress  loading  oto-  Fig.  3.19  shows  similar  results 
for  the  case  of  the  fastener  loading  on  the  hole,  while  Fig.  3.20  shows  the  combined  effect  of  far-field  loading 
as  well  as  the  fastener  loading.  As  in  the  case  of  the  central  crack  discussed  in  Sec. ,  the  stiffer  the  adhesive 
layer  and  the  stiffer  the  patch,  the  more  reduction  is  achieved  in  the  ^-factor  near  the  patched  crack-tip. 
On  the  other  hand,  a  stiff  adhesive  layer  has  a  deleterious  effect  on  the  integrity  of  the  patch  itself,  as  the 
adhesive  shear  stresses  also  increase  as  the  adhesive  stiffness  increases. 
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Figure  3.17:  Variation  of  nomalized  A:-factor  for  a  given  crack  length,  as  a  function  of  the  hole  radius(repaired 
case) 
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Figure  3.18:  Variation  of  nomalized  A:-factor,  K/b,  for  a  given  crack  length,  as  a  function  of  the  patch  stiffness 
and  adhesive  flexibility 
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Figure  3.19:  Variation  of  nomalized  /:-factor,  K/c,  for  a  given  crack  length,  as  a  function  of  the  patch  stiffness 
and  adhesive  flexibility 


Figure  3.20:  Variation  of  nomalized  k-factor,  Kja,  for  a  given  crack  length,  as  a  function  of  the  patch  stiffness 
and  adhesive  flexibility 
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Figure  3.21:  Uniaxial  fuselage  lap  joint  specimen 


Multiple-site-damage  in  a  fuselage  lap-joint 

A  model  of  a  typical  lap  joint,  and  a  typical  multiple  site  damage  near  a  row  of  fastener  holes,  are  illustrated 
in  Fig.  3.21  and  b,  respectively.  Taking  into  account  the  fastener  flexibility,  the  fastener  reaction  forces  on 
the  upper  (or  lower)  skin  can  be  determined  for  the  lap-splice  joint  configuration  of  Fig.  3.21,  through  the 
finite  element  method[Park,  Ogiso,  and  Atluri  (1992)].  Once  these  reaction  forces,  treated  as  concentrated 
forces,  are  determined,  one  may  use  the  known  elasticity  solution,  to  approximate  the  stress  field  on  the 
hole-surface  in  the  skin,  that  is  equivalent  to  these  concentrated  forces.  Under  the  action  of  these  fastener- 
interaction  stress  fields  and  the  far-field  hoop  tension,  the  stress-intensity  factors  for  multiple  cracks  near 
fastener  holes  can  be  determined  using  the  FEAM  described  in  this  section.  The  corresponding  reduction  in 
Mactors  due  to  patching  can  be  obtained  again,  using  the  FEAM,  as  described  in  this  section. 

3.1.3  Three-dimensional  analysis  of  surface-flaws  with  and  without  repairs 

In  the  FEAM  applied  to  three-dimensional  problems  of  embedded  and/or  surface  flaws  (of  elliptical  or  part- 
elliptical  shapes)  in  structural  components,  the  key  ingredient  is  the  analytical  solution  for  an  embedded 
elliptical  flaw  in  an  isotropic  elastic  solid,  the  crack-faces  being  subjected  to  arbitrary  tractions  [Atluri 
(1986);  Nishioka  and  Atluri  (1983);  Atluri  and  Tong  (1991)].  Here  we  describe  some  applications  of  the 
three-dimensional  FEAM  to  the  analysis  of  surface-flaws  in  aircraft  structural  components,  with  and  without 
composite-patch  repairs. 


Comer  flaw  near  a  counter-sunk  rivet  hole 

The  problem  of  a  counter-sunk  rivet  hole  with  a  surface-flaw  emanating  from  the  comer  is  shown  in 
Fig.  3.23.  This  is  a  typical  situation  in  an  aircraft-fuselage  lap-joint.  In  the  early  stages  of  fatigue,  this 
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2024 -T3  plate  Tolerance  =  0.025  mm 

/  =  0.040” 

’Crack  angle  =  0" 


Figure  3.22:  details  of  spark  eroded  crack  starters 


Configurations:  oH  =  0.2,  0.5, 0.8 
nt  =  2.0, 1.0,  0.5 
d//  =  0  &b/t>0.0 


Loading :  Remote  tension 
Remote  bending 


Figure  3.23:  Schematic  of  a  counter  sunk  rivet  hole  with  a  surface  crack 
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Counter  -  sunk  Rivet 

Elements:  1,430  20-node  elements 
Nodes:  6,82B 

D.O.F. :  20,484 

(Assumes  symmetry  of  plate  8<  crocks) 


90  Elements  only 
1 20  nodes ) 

202  Nodes 
606  D.O.F. 


Present  model 
(Finite  Element  Alternating  Metltodl 


Figure  3.24:  Typical  finite  element  models 


surface  flaw  grows  into  a  through-the-thickness  crack  emanating  from  the  fastener  hole  leading  ultimately 
to  the  development  of  multiple-site-damage  near  a  row  of  such  fastener  holes. 

Fig.  3.24  shows  a  comparison  of  two  different  finite-element  models:  (i)  Fig.  3.24a  shows  a  finite  element 
model  of  ( 1  /4)  of  the  plate  when  two  surface  cracks  are  assumed  to  emanate  symmetrically  from  either  side 
of  the  hole.  In  this  finite  element  model,  crack-tip  elements  are  used  in  an  explicit  numerical  modeling  of 
the  crack-tip  singularity.  Because  of  this,  this  model  involves  1430  twenty-node  finite  elements,  with  a  total 
of  20.484  degrees  of  freedom;  (ii)  Fig.  3.24b  shows  the  finite  element  model  in  the  present  FEAM.  In  this 
model,  only  one  surface  crack  is  assumed  to  emanate  from  one  side  of  the  hole;  and  the  crack-tip  singularity 
is  treated  analytically  through  the  solutions  given  in  Vijayakumar  and  Atluri  (1981)  and  Nishiokaand  Atluri 
( 1 983).  Because  of  this,  the  finite  element  model  involves  only  90  twenty-node  elements  with  a  total  of  only 
606  degrees  of  freedom.  Thus,  the  savings  in  computational  time,  as  well  as  in  the  data  preparation  time,  in 
the  present  FEAM  are  tmly  significant. 

In  Fig.  3.25,  a  comparison  of  the  results  from  a  conventional  FEM,  and  the  present  FEAM,  is  shown  for 
the  normalized  /:-factor  at  the  tip  of  the  minor  axis  of  the  quarter-elliptical  flaw,  as  a  function  of  the  ratio  of 
the  depth  of  the  hole  (without  counter-sink)  to  the  depth  of  the  plate.  The  FEAM  with  only  1575  d.o.f.  gives 
as  good  a  set  of  results  as  the  conventional  FEM  with  20,484  d.o.f. 

Three-dimensional  problem  of  a  plate  with  a  surface  crack:  effect  of  patching 

Consider  the  problem  of  a  base  plate  with  a  surface  crack,  with  a  composite  patch,  as  shown  in  Fig.  3.26. 
We  assume  that  the  sizes  of  the  base  plate  and  patch  are  infinite  (this  assumption  is  reasonable  as  long  as  the 
widths  and  lengths  of  the  plate  as  well  as  the  patch  are  about  five  times  the  length  of  the  semi-major  axis  of 
the  elliptical  surface  flaw,  as  in  Fig.  3.26),  but  that  the  thicknesses  of  the  base  plate  as  well  as  the  patch  are 
finite.  Let  the  thicknesses  of  the  base  plate  and  the  patch  be  hg  and  hp  respectively.  The  loading  is  assumed 
to  be  far-field  tension,  as  in  Fig.  3.26. 


Figure  3.25:  Variation  of  ^-factor  at  minor  axis  of  ellipse  as  a  function  of  hole  depth 


Figure  3.26:  Schematic  of  a  patched,  surface-flawed  thick  plate 
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Figure  3.27:  Linear  superposition  of  (a)  a  patched  plate  without  a  crack  subjected  to  far-field  stress,  and  (b) 
patched  plate  with  a  crack  subjected  to  crack-face  pressure 


Figure  3.28:  Linear  superposition  of  basic  problems 


The  problem  in  Fig.  3.26  can  be  solved  by  using  the  superposition  principle  as  town  in  Fig.  3.27.  In  this 
case,  the  displacement  fields  of  the  basic  problems  as  shown  in  Fig.  3.27  need  to  be  solved  for. 

Fig.  3.28a  is  the  case  wherein  uniform  stress  Oq  is  applied  on  the  crack-surface.  The  displacement  field 
for  this  problem  is  generated  by  using  an  analytical  alternating  technique,  wherein  the  analytical  solution 
for  an  infinite  body  containing  an  embedded  elliptical  flaw,  subjected  to  arbitrary  crack-face  tractions,  is 
employed.  This  analytical  solution  has  been  previously  obtained  by  Nishiokaand  Atluri  (1983). 

Fig.  3.28b  shows  the  problem  wherein  point  loads  X,Z  are  applied  symmetrically  at  the  points  (xo,zo). 
(xo,  -zo)>  (--^OiZo).  (-^0,  ~zo)  on  the  plate  with  a  semi-elliptical  surface  crack.  Let  the  displacements  at  a 
point  (A:,z)on  the  front  surface  (y  =  0)  be  expressed  as 


u  =  H\i{x,z;xo,zq)X  +  Hi2{x,z;xo,zo)Z 
W  =  H2\{x,Z\Xo,Zo)X  +  H22{x,z;Xo,Zo)Z 


The  problem  of  Fig.  3.28c  is  one  wherein  point  loads  X,Z  are  applied  symmetrically  on  an  infinite 
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Figure  3.29:  Reduction  in  ^-factor  due  to  patching  as  a  function  of  crack-depth  ratio 


orthotropic  plate.  This  solution  has  been  reported  on  previously.  Let  the  displacements  at  a  point  (jc,z)  for 
this  problem  be  expressed  as: 


=  A:ii  (x,z;.ro,zo)^  +  /^i2(v.r;  V:i.-o)Z 
=  A'21  ^  +  ^^22  (  v. '-2:0)2 

Let  the  displacements  U5  and  w5  denote  the  total  displacements  from  problems  3.28A  and  B.  The 
compatibility  conditions  between  the  base  plate  and  the  composite  patch  are  expressed  as 

ir'  -  vv'''  -  wP 

=  j.  •  L 

where  F  is  the  adhesive  flexibility,  defined  as  earlier,  as 

f  = 

~  Ga  "  8G,, 

By  using  these  compatibility  conditions,  integral  equations  are  derived  for  Xxz  and  lyz-  Once  these  are 
solved  for,  the  /:-factors  are  obtained  in  a  manner  described  in  earlier  sections. 

Fig.  3.29  shows  the  effect  of  the  thickness  of  the  base  plate  on  the  normalized  stress-intensity  factors 
for  both  the  unrepaired  and  repaired  cases.  When  the  base  plate  thickness  is  small,  about  50%  reduction 
in  k  -factor  is  obtained.  As  the  base  plate  becomes  thicker,  however,  the  effect  of  the  patch  is  reduced  (to 
about  40%),  and  the  A:-factor  converges  to  a  constant  value.  (Note  that  in  the  case  of  the  surface-flaw,  the 
jt-factor  that  is  considered  is  at  the  intersection  of  the  crack-front  with  the  front  face  of  the  base  plate).  From 
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Figure  3.30:  Effects  of  patch  stiffness  and  adhesive  flexibility  on  the  /:-factor  of  a  surface  flaw 


Fig.  3.29  it  is  also  noted  that  if  one  uses  a  very  stiff  patch,  a  better  reduction  in  the  yt-factor  is  obtained. 
However,  a  stiff  patch  may  increase  the  adhesive  shear  stresses. 

Fig.  3.30  shows  the  effect  of  patch-stiffness  parameters  and  adhesive  flexibility  parameters  on  the  nor¬ 
malized  /:-factors. 

Fig.  3.31  shows  the  effect  of  patch  stiffness  and  adhesive  flexibility  on  the  adhesive  shear  stresses,tyz/ao 
at  A  =  0,  y  =  0  and  z  =  0  .  Finally  Fig.  3.32  shows  the  effect  of  the  patch-stiffness  and  adhesive  flexibility 
on  the  values  of 

Repair  of  a  surface  crack  near  a  hole 

The  problem  is  shown  schematically  in  Fig.  3.33.  To  repair  the  crack,  a  patch  of  width  2Wp  and  height  IHp 
is  applied,  such  that  the  center  of  the  patch  coincides  with  the  center  of  the  hole. 

The  basic  component  problems  to  be  analyzed  are  shown  in  Fig.  3.34a-c.  For  problem  3.34a,  the  dis¬ 
placement  fields  on  the  front  surface  and  the  stress  intensity  factors  can  be  obtained  by  using  the  3D  finite 
element  alternating  technique  [Atluri  (1986)]. 

As  for  problem  3.34b,  when  a  point  force  (X,Z)  is  applied  at  the  point  (ao,zo)  and  a  symmetrical  point 
load  (X,  -Z)  is  applied  at  (xq,  -zo),  once  again  the  3D  finite  element  alternating  technique  can  be  used. 

In  problem  3.34c,  the  displacement  field  in  the  patch  is  determined  as  before. 

The  finite  element  mesh  used  for  problems  3.34a  and  b  are  shown  in  Fig.  3.35.  The  finite  element  mesh 
used  for  problem  3.34c  is  shown  in  Fig.  3.36. 
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Figure  3.32:  Effects  of  patch  stiffness  and  adhesive  flexibility  on  the  tensile  stress  in  the  patch 


Figure  3.33:  Schematic  of  a  surface  flaw  near  a  fastener  hole 


Figure  3.36:  FEM  mesh  for  patch  plate,  64  2D  8  node  isoparametric  elements,  225  nodes 
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Figure  3.37:  Variation  of  S.I.F  near  the  border  of  a  surface  flaw  near  a  fastener  hole,  with  and  without 
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The  variations  of  the  stress  intensity  factor  along  the  crack  front,  with  and  without  repair  patch,  are 
shown  in  Fig.  3.37.  Once  again,  the  significant  reduction  in  the  crack-growth  retardation,  achieved  by  the 
application  of  the  patch,  is  evident. 

Simplified  computational  procedures  for  the  analysis  of  structural  integrity  problems  in  aging  airplanes, 
with  and  without  composite-patch  repairs,  have  been  developed,  and  their  use  is  illustrated. 

3.2  Patch  repair  and  curved  crack 

3.3  Analysis  Of  Repaired  Cracks  in  Pressurized  Aircraft  Fuselages 

Repairs  to  pressurized  fuselages  have  been  done  traditionally  by  using  mechanical  doublers.  While  bond¬ 
ing  of  metal  structures  and  bonded  repair  techniques  have  existed  for  about  50  years  now,  bonded  repairs 
have  not  gained  as  much  acceptance  as  the  mechanical  doubler  repairs.  This  is  mainly  due  to  the  disap¬ 
pointing  early  experience  with  bonding  and  the  lack  of  widespread  understanding  of  bonded  repairs  among 
technicians.  This  provides  the  motivation  for  the  development  of  effective  analysis  techniques  that  can  be 
implemented  on  low  end  workstations.  In  this  investigation,  the  mechanical  doubler  repairs  have  been  com¬ 
pared  with  bonded  repairs  in  the  repair  of  cracks  in  the  fuselage  skin.  To  analyze  the  effect  of  global  loading 
(viz.  pressure  loading  in  the  fuselage  shell),  on  the  essentially  local  feature  of  repair,  a  hierarchical  approach 
has  been  used.  This  allows  the  fuselage  to  be  modeled  with  increasing  detail  over  smaller  regions.  At  the 
global  level,  the  fuselage  is  modeled  as  a  shell,  with  stiffeners  modeled  as  beams,  and  with  flexible  fasteners. 
In  the  intermediate  stage  the  degree  of  modeling  is  improved.  A  section  of  the  fuselage  from  the  global  stage 
is  now  analyzed,  with  the  stiffeners  being  modeled  with  shell  elements.  The  kinematic  boundary  conditions 
applied  to  this  smaller  intermediate  region  are  obtained  from  the  global  solution.  The  stress  results  from  this 
intermediate  zone  are  later  applied  to  a  still  smaller  zone  that  is  analyzed  in  the  local  stage;  and  this  local 
analysis  is  used  to  evaluate  the  stress  intensity  factors  for  repaired  or  unrepaired  cracks.  The  local  stage 
uses  a  2D  finite  element  alternating  method  (FEAM)  to  model  the  crack.  At  this  stage,  the  bending  effects 
caused  by  a  bonded  patch  on  one  side  of  the  skin  are  ignored.  TIt:  assumption  is  considered  sufficient 
as  the  substructure  of  the  fuselage  will  restrict  the  effect  of  bending  trom  being  of  any  significance.  The 
local  analysis  accounts  for  the  non-linearity  of  the  adhesive  material  by  modeling  it  as  an  elastic  perfectly 
plastic  material.  A  large  difference  in  the  extensional  stiffness  of  patch  and  the  skin  results  in  large  strains 
being  induced  in  the  adhesive.  Therefore,  the  necessity  to  model  the  adhesive  as  an  elastic  perfectly  plas¬ 
tic  material  becomes  more  critical  when  the  patch  extensional  stiffness  is  large  when  compared  to  the  skin 
extensional  stiffness.  Using  this  approach  the  stress  redistribution  caused  by  mechanical  doubler  repair  is 
compared  with  that  caused  by  composite  patch  repairs.  The  residual  strength  and  fatigue  life  of  fuselages 
with  composite  patch  repairs  have  been  compared  with  those  of  the  un-repaired  case.  The  stress  intensity 
factors  have  been  evaluated  by  considering  the  adhesive  as  elastic  material  and  as  an  elastic  perfectly  plastic 
material.  Parametric  studies  of  composite  repairs  are  conducted  using  this  methodology. 

Most  of  the  present  day  repairs  to  cracked  aircraft  fuselages  are  performed  by  using  mechanical  doublers. 
The  cracked  portion  is  cutout,  a  sheet  is  placed  in  the  region  thus  created;  and  finally,  another  sheet  is  riveted 
to  this  cutout  portion  and  the  skin. 

From  the  Second  World  War,  successful  applications  of  adhesive  bonding  have  been  made  in  aircraft 
structures.  The  Fokker  F-27  aircraft  is  an  example  of  successfully  operating  aircraft  that  has  employed 
adhesive  bonding  for  primary  structure  [Baker  and  Jones  ( 1 988)].  While  such  success  stories  exist ,  doubts 
about  adhesive  bonded  repairs  still  exist.  In  this  section  we  study  such  bonded  repairs  and  compare  them 
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with  mechanical  doubler  repairs. 

There  are  several  methods  for  the  analysis  of  patched  cracks.  They  can  be  broadly  divided  as 

1.  Analytical 

2.  Finite  Element  Approach 

The  analytical  approach  of  Rose  (1981)  is  an  elegant  method  based  on  Hart-Smith  (1974)  theory  of 
bonds,  elastic  inclusion  analogy,  and  on  some  simplifying  assumptions.  Fredell  (1994),  in  his  exhaustive 
thesis,  has  extended  this  analysis  to  include  thermal  effects.  He  has  also  carried  out  an  evaluation  of  mechan¬ 
ical  doubler  repairs.  Erdogan  and  Arin  (1972)  have  used  integral  equations  approach  to  study  the  patched 
cracks.  The  assumptions  of  Erdogan  and  Arin  were  subsequently  used  by  Ko  (1978)  and  Hong  and  Jeng 
(1985)  in  the  analysis  of  sandwich  plates  with  part-through  crack. 

Jones  and  co-workers  [Jones  and  Callinan  (1979)],  Mitchell,  Wooley,  and  Chivirut  (1975);  Chu  and  Ko 
(1989)  used  finite  element  method  to  study  patched  cracks.  Park,  Ogiso,  and  Atluri  (1992)  have  used  an 
integral  approach  combined  with  FEAM  to  estimate  the  stress  intensity  factors  for  patched  panels.  Tam  and 
Shek  (1991)  have  combined  boundary  element  method  approach  (for  the  plate)  and  finite  element  method 
(for  the  patch)  to  estimate  the  stress  intensity  factors.  Other  work  in  this  area  includes  Atluri  and  Kathiresan 
(1978),  Sethuraman  and  Mathi  (1989),  and  Kan  and  Ratwani  (1981). 

In  most  of  these  approaches,  only  patches  of  infinite  size,  or  very  narrow  strip  type  patches,  or  infinite 
sheet  cases  are  considered.  All  of  these  cases  are  valid  only  for  flat  sheets.  The  loading  for  all  these  analyzes 
are  hoop  stresses  evaluated  from  basic  thin  shell  theory.  While,  in  most  cases  this  is  a  good  approximation, 
this  does  not  take  into  account  the  stress  re-distributions  due  to  the  curvature  and  due  to  the  presence  of 
stiffeners. 

The  present  work  addresses  the  problem  of  the  study  of  patched  cracks  in  actual  commercial  airliner 
fuselages.  A  hierarchical  modeling  strategy  is  used  to  study  the  repairs  to  fuselages.  The  methodology 
is  as  follows.  The  whole  of  the  fuselage  is  studied  in  the  first  stage.  In  the  next  stage,  a  portion  of  this 
is  studied  with  an  increase  in  complexity  of  the  model.  Finally  a  small  portion  around  the  crack  and  the 
repairs  is  studied  to  estimate  the  stress  intensity  factors.  The  procedure  used  to  analyze  the  local  zone 
is  the  finite  element  alternating  method.  This  methodology  to  study  repairs  to  fuselages.  The  stress  re¬ 
distributions  caused  by  mechanical  doubler  repairs  is  contrasted  with  those  by  patch  repairs.  The  effect 
of  various  parameters  such  as  the  patch  dimensions,  patch  material,  patch  thickness,  adhesive  material  and 
adhesive  thickness  on  patch  efficiency  is  studied. 

3.3.1  Methodology 

The  problem  of  cracks  and  their  repairs  is  a  localized  phenomenon.  The  known  loading  conditions  on 
the  fuselage  are  essentially  global  -  viz.  the  pressure  loading.  Therefore  the  analysis  requirements  are 
conflicting  -  a  global  loading  and  an  essentially  local  phenomenon.  This  necessitates  a  hierarchical  modeling 
strategy  (Starnes  and  Britt  1991).  Fig.  3.3.1.  This  allows  an  increasing  level  of  complexity  and  fidelity  in 
the  modeling  without  a  prohibitively  high  computational  effort. 

The  global  model  is  a  linear  elastic  shell  finite  element  analysis.  In  this  analysis  the  stiffeners  are 
modeled  as  beams.  The  fuselage  is  modeled  with  shell  elements.  The  fasteners  are  modeled  as  shear 
springs.  In  the  analysis  of  the  cracked  unrepaired  fuselage,  the  crack  is  modeled  in  the  global  analysis  only 


Global  Analysis 


Local  Analysis 


Figure  3.38;  Schematic  of  the  Hierarchical  Modeling  Strategy 


38 


crudely,  i.e.,  via  unconnected  finite  elements  while  ignoring  the  crack-tip  singular  fields.  In  the  analysis  of 
a  composite  patch-repaired  fuselage,  it  is  assumed  that  the  bending  caused  by  the  patch  dies  down  and  at 
this  global  stage  of  analysis,  it  can  be  considered  as  an  uncracked  fuselage.  Therefore,  the  repaired  fuselage 
is  considered  as  an  intact  fuselage  in  the  global  level  of  analysis.  In  both  cases,  symmetry  conditions 
approximating  a  full  fuselage  are  applied  to  the  region  analyzed.  As  the  shell  length  considered  is  small 
when  compared  to  the  length  of  actual  fuselages,  the  axial  ends  are  restrained  from  motion  in  the  axial 
direction. 

A  portion  of  the  fuselage  in  the  circumferential  as  well  as  longitudinal  directions  is  now  taken  and 
analyzed  in  the  intermediate  model.  In  this  model,  the  fidelity  is  increased  by  modeling  the  stiffeners  with 
shell  elements.  The  fasteners  are  modeled  by  short  beams.  The  essential  boundary  conditions  for  this  model 
are  obtained  from  the  global  analysis.  In  both  of  the  above  models  the  loading  is  internal  pressure. 

A  smaller  portion  of  the  fuselage  that  contains  the  crack  and  repair  is  then  analyzed  in  the  local  analysis. 
As  the  region  is  small,  the  pressure  loading  and  the  curvature  are  disregarded.  The  loading  on  this  region 
are  the  tractions  obtained  from  the  intermediate  model.  The  local  analysis  proceeds  in  two  stages.  In  the 
first  stage,  the  crack  is  explicitly  modeled  with  ordinary  finite  elements,  (i.e  ,  without  paying  attention  to 
the  singular  fields)  and  the  stresses  exerted  by  the  composite  patch  on  the  metallic  sheet  are  deduced.  The 
mesh  for  this  analysis  is  not  detailed  as  the  aim  at  this  stage  is  to  obtain  the  stresses  exerted  by  the  patch  and 
adhesive  on  the  sheet  and  not  to  model  the  crack  tip  singularity.  Once  the  patch  stresses  are  deduced,  the 
finite  element  alternating  method  is  used  to  obtain  the  stress  intensity  factors. 


Global  Analysis 


The  stiffeners  at  this  stage  are  modeled  as  beams.  The  beam  element  is  obtained  by  collapsing  the  shell 
element  described  above.  As  the  fastener  flexibility  is  crucial  (Swift  (1974),  (1984)),  the  fasteners  are 
modeled  by  flexible  shear  springs.  The  stiffness  for  this  fastener  element  is  given  as  in  Swift  (1974,  1984) 
by 


K  = 


(3.58) 


where  Es  is  the  sheet  Young' s  Modulus,  D  is  the  rivet  diameter,  t\  and  t2  thicknesses  of  joined  sheets, 
4  =  5.0  for  A1  rivets  and  1 .66  for  steel  fasteners,  and  C  =  0.8  for  A1  rivets  and  0.86  for  steel  fasteners.  The 
rivet  holes  are  not  explicitly  modeled. 

The  crack  is  modeled  with  split  nodes.  There  is  no  explicit  modeling  of  the  crack  tip  singularity. 


Intermediate  Analysis 

The  intermediate  analysis  is  also  a  linear  finite  element  analysis  and  is  also  based  on  the  aforementioned 
principles.  The  element  used  at  this  stage  is  a  four-noded  six  degree  of  freedom  per  node  element  (Rankin 
and  Brogan(1991)).  In  the  analysis  of  shells  a  two  parameter  specification  of  the  the  director  field  about 
the  mid-surface  is  sufficient  if  the  mid-surface  is  smooth.  In  the  presence  of  shell  intersections,  this  is 
not  sufficient.  A  standard  way  of  handling  such  intersections  is  the  usage  of  shell  elements  with  drilling 
degrees  of  freedom  (  See  Simo  ( 1 993)  and  Fox  and  Simo  ( 1 992)  for  example).  The  element  used  is  one  such 
element.  Therefore,  the  facets  that  arise  as  a  result  of  modeling  the  stiffeners  with  shells  elements  can  be 
handled.  This  element  uses  an  incompatible  displacement  field  with  a  cubic  variation  of  the  bending  field 
and  a  linear/cubic  variation  of  the  in-plane  field. 
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In  this  model,  the  stiffeners  are  also  modeled  with  the  same  shell  elements.  The  fasteners  are  modeled 
as  short  beams.  As  in  the  global  model  the  rivet  holes  are  not  explicitly  modeled.  Again  the  cracks  are 
modeled  with  split  nodes  and  there  is  no  explicit  modeling  of  the  stress  singularity. 


Local  Analysis 

The  local  analysis  is  a  two  stage  analysis.  They  are 

1.  Evaluation  of  the  stresses  exerted  by  the  adhesive  and  patch  on  the  sheet,  using  a  coarse  mesh. 

2.  The  stresses  obtained  from  stage  1  are  applied  as  body  forces  on  the  base  sheet,  and  the  finite  element 
alternating  method  is  used  to  find  the  stress  intensity  factor. 


In  step  1,  a  traditional  finite  element  methodology  is  used  to  deduce  the  stresses  exerted  by  the  patch. 
As  the  crack  tip  is  not  meshed  for  the  singularity,  a  coarse  mesh  is  sufficient  for  this  purpose.  The  sheet  and 
the  patch  are  modeled  with  eight  noded  2D  elements.  The  adhesive  is  modeled  by  elements  as  in  Jones  and 
co-workers  [Jones  and  Callinan  (1979)],  Mitchell,  Wooley,  and  Chivirut  (1975).  This  element  is  obtained 
by  assuming  a  constant  shear  stress  in  the  adhesive.  Under  the  present  conditions  the  shear  stress  can  be 
obtained  from 


'^XZ  — 


(3.59) 


where  u  and  v  are  the  displacements  in  the  x  direction  and  y  direction  respectively,  t,  g  are  the  thickness  and 
shear  modulus  respectively.  The  subscripts  a,  s  and  p  refer  to  the  adhesive,  sheet  and  the  patch  respectively. 

Large  strains  are  induced  in  the  adhesive  when  the  difference  between  the  patch  stiffness  and  sheet  stiff¬ 
ness  is  large.  Therefore,  this  results  in  the  yielding  of  the  adhesive.  This  reduces  the  patched  efficiency,  and 
a  linear  elastic  analysis  would  grossly  overestimate  the  efficiency  and  would  be  anti-conservative.  There¬ 
fore,  the  adhesive  is  modeled  as  an  elastic  perfectly  plastic  material .  \s  there  is  no  unloading  ( crack  growth 
is  not  considered),  a  deformation  theory  of  plasticity  is  considered  sufficient.  An  initial  stress  algorithm 
[Nayak  and  Zienkiewicz  (1972)]  has  been  implemented.  In  view  of  the  model  as  in  Jones  and  co-workers 
[Jones  and  Callinan  (1979)],  Mitchell,  Wooley,  and  Chivirut  ( 1 975),  the  Von  Mises  yield  criterion  reduces 
to 


(3.60) 


The  algorithm  for  this  is  as  follows 


1.  A  load  step  is  taken  and  the  displacements  are  found. 

2.  The  stresses  are  computed  at  each  Gauss  Point  and  yielding  is  checked  for. 

3.  The  residual  load  is  computed. 

4.  The  displacements  for  this  residual  load  are  then  computed. 

5.  Repeat  steps  2,  3  and  4  for  convergence.  Convergence  is  reached  when  there  is  no  yielding  at  any 
Gauss  point. 
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6.  Repeat  the  steps  1,2, 3,4  and  5  for  more  load  steps  till  the  desired  load  is  reached. 

In  finite  element  alternating  method,  the  stresses  in  the  uncracked  body  are  first  analyzed,  by  a  traditional 
finite  element  method,  for  the  given  system  of  external  loading.  To  model  the  crack,  the  tractions  at  the 
locations  of  the  crack  in  an  otherwise  uncracked  body  must  be  erased. 

Software 

The  above  methodology  has  been  implemented  in  a  software  based  on  several  in-house  and  commercial 
finite  element  codes.  The  global  model  has  been  implemented  in  an  in-house  software  called  SOFRAC 
[Singh,  Park,  and  Atluri  (1994).  The  results  from  this  analysis  are  taken  and  applied  to  the  intermediate 
analysis.  At  this  point  it  is  relevant  to  point  out  the  methodology  adopted  to  interpolate  the  displacements 
between  the  models.  While  the  fuselage  shell  has  been  modeled  with  shell  elements  in  both  the  models,  the 
stiffeners  in  the  global  model  was  modeled  as  beams;  but  are  modeled  as  shells  in  the  intermediate  analysis. 
Therefore,  a  direct  application  of  the  displacements  at  the  relevant  locations  becomes  impossible.  Instead, 
using  the  rotations  and  displacements  in  the  beams  that  constitute  the  stiffeners  in  the  global  analysis,  the 
displacements  at  the  nodes  of  the  shell  elements  that  constitute  the  stiffeners  in  the  intermediate  model  are 
obtained  by  extrapolation.  It  is  noted  at  this  point  this  could  lead  to  a  stress  mismatch  in  the  boundary  of  the 
intermediate  model  as  the  idealization  is  different. 

The  intermediate  model,  with  the  kinematic  boundary  conditions  obtained  as  illustrated  above,  is  now 
solved  by  using  STAGS  [Almroth,  Brogan,  and  Stanley  ( 1 986)].  The  pre-processor  used  to  generate  the  input 
files  for  STAGS  analysis  is  an  in-house  pre-processor  that  carries  out  the  aforementioned  extrapolations. 

The  solution  from  the  above  is  now  used  as  the  applied  boundary  condition  for  the  local  analysis  wherein 
only  a  very  small  portion  of  the  skin(  with  the  patch)  are  modeled.  As  indicated  above,  the  local  analysis 
is  based  on  2-D  finite  element  alternating  method.  Therefore,  only  the  in-plane  stresses  are  used.  The  out- 
of-plane  stresses  are  neglected.  As  the  substructure  of  the  structure  would  restrict  any  severe  bending  this 
is  believed  to  be  a  good  approximation.  The  local  zone  is  implemented  in  an  in-house  software  PAN-MD. 
The  input  files  for  this  analysis  are  generated  by  an  in-house  post-processor  for  STAGS.  The  entire  analysis 
from  the  global  to  local  has  been  automated  in  a  PERL  script. 

3.3.2  Analysis  of  Repaired  Cracks 

The  fuselage  as  described  by  Tab.  3.3.2  is  used  for  the  various  test  cases  studied.  The  stringers  used  are  of 
hat  cross  section  and  the  frames  are  of  Zee  cross  section.  The  stiffener  properties  are  given  in  Tab.  3.3.2. 
It  is  noted  that  the  frame  is  attached  to  the  sheet  on  the  tearstraps.  This  is  considered  as  an  approximation 
to  attaching  the  frame  to  the  sheet  with  shear  clips.  In  view  of  the  unavailability  of  the  exact  stiffener 
dimensions  this  is  considered  sufficient.  The  crack  is  assumed  to  be  a  mid-bay  crack.  The  rivets  used  are 
NAS1097DD6  rivets  of  0.15625  in  diameter. 

The  dimensions  of  the  different  models  are  given  in  Tab.  3.3.  The  finite  element  meshes  used  for  various 
models  are  shown  in  Fig.  3.39,  Fig.  3.40  and  Fig.  3.41 . 

Patch  Repairs 

The  above  fuselage  is  repaired  with  a  patch  of  dimensions  5  in  X  3  in  as  shown  in  Fig.  3.42.  The  patch  is 
0.04  in  thick  and  is  made  of  Boron  Epoxy.  The  adhesive  is  0.004  in  thick  and  has  a  shear  modulus  of  104.73 


Tab: 


e  3.1:  Properties  of  The  Fuselage  Shell 


Internal  Pressure 

9.0  psi 

Radius 

78  in 

Thickness 

0.036  in 

Material 

A1  2024T3 

Young' s  Modulus 

10,500  ksi 

Poisson’s  Ratio 

0.33 

Yield  Point 

47.0  ksi 

Table  3.2:  Stiffener  Properties 


Rivet  Specification 

NAS1097DD6 

Rivet  Diameter 

0.15625  in 

Stringer  Spacing 

9.25  in 

Stringer  Area 

0.0384  in^ 

Stringer  Moment  of  Inertia 

0.162  in'* 

Axis  to  skin  distance 

0.465  in 

Frame  Spacing 

20.0  in 

Frame  Area 

0.1248  in- 

Frame  Moment  of  Inertia 

0.3271  m" 

Axis  to  skin  distance 

1.148  in 

Tear  Strap  Thickness 

0.036  in 

Tear  Strap  Width 

2  in 

Table  3.3:  Dimension  of  different  fuselage  models 


Global  Region  Length 

100  in 

Global  Region  Arc  Length 

46.25  in 

Intermediate  Region  Length 

40  in 

Intermediate  Region  Arc  Length 

18.5  in 

Local  Region  Length 

10  in 

Local  Region  Width  (Arc  Length) 

8  in 
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Figure  3.39:  Mesh  for  Global  Analysis 


Figure  3.40:  Mesh  for  Intermediate  Analysis 
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Figure  3.41:  Mesh  for  Local  Analysis 


ksi  with  5.250  ksi  being  the  Yield  Point  in  shear. 

To  analyze  the  effect  of  the  patch,  the  crack  length  was  varied  from  1  in  to  16  in.  The  variation  of  the 
stress  intensity  factor,  K[,  with  the  crack  length  is  shown  in  Fig.  3.43.(Unless  otherwise  noted  all  the  stress 
intensity  factors  are  in  ksiy/Jn  and  the  crack  lengths  are  in  in.)  The  efficiency  of  crack  patching  is  plotted  in 
Fig.  3.44.  The  efficiency  is  defined  as 

^^^unpatched  ~  ^patched  , 

= - ^ - X  100 

^  unpatched 

The  asymptotic  nature  of  the  stress  intensity  factor,  as  in  Fredell  (1994)  and  Baker  and  Jones  (1988)  is 
clearly  observed.  Unlike  the  methodology  of  Rose  (1981)  and  Fredell  (1994),  the  present  methodology  can 
analyze  cracks  longer  than  the  width  of  the  patch.  It  is  seen  that  while  the  stress  intensity  factors  show  a 
significant  increase  as  the  crack  emerges  from  the  patch,  it  still  demonstrates  an  asymptotic  behavior.  In 
these  two  situations,  viz.  crack  shorter  than  the  width  of  the  patch  and  longer  than  the  width  of  the  patch, 
the  patch  transfers  load  in  significantly  different  ways.  When  the  crack  is  shorter  than  the  width  of  the 
patch,  most  of  the  load  is  transfered  by  the  edges  of  the  patch  parallel  to  the  crack.  This  is  seen  from  the 
Fig.  3.46,  where  the  effective  stress  in  the  adhesive  has  been  plotted.  As  the  crack  grows  longer  than  the 
width  of  the  patch,  the  regions  near  the  crack  tips  also  start  to  transfer  the  load  and  the  patch  transfers  the 
load  as  in  bonded  overlap  joints.  This  is  shown  in  the  effective  stress  in  the  adhesive  contour  plot  Fig.  3.47. 
This  increase  in  the  effective  stress  lead  to  a  final  yielding  of  the  adhesive.  Once  the  adhesive  yields,  a  very 
significant  reduction  in  the  effectiveness  in  seen.  (Fig.  3.44  and  Fig.  3.45). 

Patch  Repairs  and  Mechanical  Doubler  Repairs 

A  mechanical  doubler  repair.  Fig.  3.48,  for  the  above  mid-bay  cracked  fuselage  is  also  analyzed.  The  axial 
(Fig.  3.49)  and  hoop  stress  (Fig.  3.50)  resultants  for  the  doubler  repaired  fuselage  has  been  compared  with 


Figure  3.42:  Geometry  of  the  Local  Zone  showing  the  patch 
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Figure  3.43;  Stress  Intensity  factors  Vs.  Crack  Length 
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Figure  3.44:  Efficiency  of  Crack  Patching 


Figure  3.45:  Effect  of  Plasticity 
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Figure  3.46:  Effective  Stress  in  the  Adhesive  for  Crack  Length  3  in 
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Figure  3.47:  Effective  Stress  in  the  Adhesive  for  Crack  Length  15  in 
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Figure  3.48:  Geometry  of  the  doubler  repair 
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Figure  3.49:  Axial  Stress  Resultant  for  Doubler  Repair 


that  of  the  composite  patch  (Fig.  3.51  and  Fig.  3.52)  repaired  fuselage  and  the  intact  fuselage  (Fig.  3.53  and 
Fig.  3.54).  As  can  be  seen,  a  patch  repair,  despite  the  presence  of  the  singularity  involves  a  much  lower 
value  stress.  In  the  doubler  repair  case,  the  cutout  provides  a  site  for  stress  concentration.  Therefore,  the 
maximum  hoop  stress  in  the  doubler  repair  case  increases  from  19.1  ksi  in  the  intact  fuselage  to  30  ksi. 
Whereas,  in  the  fuselage  repaired  with  a  patch,  the  hoop  stress  went  up  marginally  to  1 9.4  ksi.  A  similar 
pattern  is  also  seen  in  the  axial  stresses  where  the  stresses  went  up  from  6.5  ksi  to  1 1 .7  ksi. 

Effect  of  Patch  Geometry 

To  study  the  effect  of  the  patch  overlap  length  of  the  patch  (the  length  along  the  y  direction),  it  was  varied 
from  1.5in  to  3.5  in.  The  results  are  in  Fig.  3.55.  The  patch  is  5  in  wide  and  0.04  in  thick  and  is  made  of 
Boron/Epoxy  composite.  The  adhesive  is  0.004  in  thick  and  has  a  shear  modulus  of  104.73  ksi. 
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Figure  3.54:  Hoop  Stress  Resultant  for  Intact  Fuselage 
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Figure  3.55:  Effect  of  Patch  Overlap  Length  on  Stress  Intensity  Factor 
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Figure  3.56:  Effect  of  Patch  Overlap  Length  on  Stress  Intensity  Factor 


When  the  crack  is  longer  than  the  width  of  the  patch,  the  patch  serves  to  transfer  the  load  across  the  crack 
and  reduce  the  ligament  stress  near  the  crack  tip.  Therefore,  a  longer  patch,  having  a  longer  overlap  length 
and  therefore  larger  area  transfers  this  stress  better  and  therefore  perform  better  than  the  shorter  patches 
(Fig.  3.55  and  Fig.  3.56).  Whereas,  when  the  crack  is  shorter  than  the  width  of  the  patch,  the  longer  patches 
perform  poorly  (Fig.  3.57).  The  composite  system  of  the  patch,  adhesive  and  the  sheet  attracts  parasitic 
load  from  other  regions  in  the  sheet,  [Baker  and  Jones  (1988).  When  the  crack  is  shorter  than  the  width  of 
the  patch,  this  parasitic  load  attaction,  increases  the  load  near  the  crack  tip.  Larger  patches  attract  larger 
parasitic  loads,  therefore  they  perform  poorly  when  compared  to  shorter  ones.  When  the  crack  is  longer, 
this  parasitic  load  attraction  does  not  increase  the  load  near  the  crack  tip,  therefore  a  similar  behavior  in  not 
seen  in  those  cases. 

The  patch  width  is  varied  from  4  in  to  8  in.  The  results  are  shown  in  Fig.  3.58.  The  patch  is  3  in  long, 
and  0.04  in  thick  and  is  made  of  Boron/Epoxy.  The  adhesive  properties  remain  the  same  as  above.  As  was 
illustrated  earlier,  the  stress  intensity  factors  shows  two  asymptotic  values,  one  obtained  when  the  crack  is 
shorter  than  the  width  of  the  patch  and  the  other  when  the  crack  is  longer  than  the  width  of  the  patch.  When 
the  crack  is  shorter  than  the  width  of  the  patch,  the  asymptotic  value  reached  does  not  vary  much  with  the 
patch  width.  While  wider  patches  cause  the  second  asmptotic  value  to  be  reached  at  a  longer  crack  length 
(the  cracks  emerge  from  the  patch  later  with  wider  patches),  the  asymptotic  value  reached  does  not  aprear  to 
change.  Therefore  the  narrower  patches  also  seem  to  perform  adequately.  Wider  patches,  because  of  larger 
areas,  have  better  damage  tolerance. 

The  patch  thickness  was  varied  from  0.01  in  thick  to  0.04  in.  The  patch  is  5  in  wide  and  3  in  long  and  is 
made  of  Boron/Epoxy.  As  expected,  thicker  patches  perform  better  (Fig.  3.59). 


Effect  of  Patch  Material 

To  evaluate  different  materials  as  possible  candidates  for  patching  materials  several  are  investigated,  follow¬ 
ing  Fredell  (1994).  The  properties  of  the  materials  are  given  in  Tab.  3.3.2.  In  each  of  the  cases  the  patch 
is  5  in  wide  and  3  in  long  and  0.04  in  thick.  The  adhesive  shear  modulus  is  104.7  ksi  and  the  thickness 
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Figure  3.57:  Eifect  of  Patch  Overlap  Length  on  Stress  Intensity  Factor 


Figure  3.58:  Effect  of  Patch  Width  on  Stress  Intensity  Factor 
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Figure  3.59:  Effect  of  Patch  Thickness  on  Stress  Intensity  Factor 


Table  3.^ 

1:  Properties  of  Various  Patching  Materials  Considered 

Material 

El 

El 

Vi2 

G\2 

Boron/Epoxy 

3.02  X  lO’  psi 

3.69  X  10^  psi 

0.1677 

1.05  X  10^  psi 

Carbon/Epoxy 

1.97  X  lO’  psi 

1.814  X  10^  psi 

0.30 

1.015x10^  psi 

A1  2024-T3 

1.05  X  10^  psi 

1.05  X  10^  psi 

0.32 

3.977  X  10^  psi 

GLARE 

9.565  X  10^  psi 

7.358  X  10^  psi 

0.33 

2.685  X  10^  psi 

considered  is  0.004  in.  The  Yield  Stength  of  the  adhesive  is  taken  as  5254  psi.  The  results  are  depicted  in 
Fig.  3.60. 

The  composite  patches  perform  the  best,  followed  by  A1 2024-T3,  and  finally  the  metal  laminate  GLARE. 
When  the  crack  is  under  the  patch,  the  difference  is  marked  .  But  once  the  crack  becomes  longer  than  the 
width  of  the  patch,  the  difference  is  marginal.  When  the  crack  is  shorter  than  the  width  of  the  patch,  the 
induced  strains  cause  larger  stress  to  be  transfered  in  stiffer  patches.  Therefore,  the  load  near  the  crack  tip  is 
reduced  and  the  stress  intensity  factor  reduction  is  more  pronounced  for  stiffer  patches.  On  the  other  hand 
when  the  crack  is  longer  than  the  width,  the  patch  transfers  the  load  as  in  a  bonded  lap  joint  and  the  load 
transferred  is  the  ligament  stress  and  has  little  relationship  with  the  stiffness.  This  is  confirmed  by  the  fact 
that  when  the  crack  is  extremely  long,  the  stress  intensity  factor  for  patched  cracks,  with  patch  of  various 
materials  is  almost  the  same. 

Effect  of  Adhesive  Properties 

With  the  effect  of  the  patch  material  and  geometry  determined  the  effect  of  the  adhesive  is  now  examined. 
With  the  boron-epoxy  patch  of  5  in  X  3  in  X  0.04  in,  the  adhesive  thickness  is  varied  from  0.001  to  0.005 
for  a  crack  of  length  8  in.  The  results  are  shown  in  Fig.  3.61.  Thin  adhesive  layers  perform  better  than  thick 
adhesive  layers  because,  a  thick  adhesive  makes  the  composite  system  of  the  sheet,  adhesive  layer  and  patch 
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Figure  3.60:  Comparison  of  Different  Patching  Materials 


softer.  Often,  a  thick  adhesive  layer  is  needed  for  good  durability  as  a  thick  adhesive  will  experience  lesser 
strains.  As  can  be  seen  from  Fig.  3.62  the  detrimental  effect  of  having  a  thick  adhesive  is  not  severe. 

The  shear  modulus  of  the  adhesive  is  varied  for  the  same  patch  dimensions  and  for  0.004  in  thick  adhe¬ 
sive.  The  result  is  shown  in  Fig.  3.62.  As  can  be  seen,  stiffer  adhesives  were  found  to  perform  better.  While, 
the  conditions  considered  were  insufficient  for  adhesive  yield,  a  stiffer  adhesive  would  result  in  the  yield 
point  being  reached  earlier.  Thus,  the  material  would  yield  and  reduce  the  patch  effectiveness  as  illustrated 
earlier.  This  being  the  case,  it  is  also  noted  that  the  advantage  in  using  a  stiffer  adhesive  (in  terms  of  stress 
intensity  reduction)  is  not  significant. 


3.3.3  Deficiencies  of  the  methodology 

Cracks  on  fasteners  rivet  rows  can  be  studied  using  a  different  local  analysis  [as  in  Park,  Ogiso,  and  Atluri 
(1992)].  As  the  methodology  is  based  on  a  finite  element  alternating  method,  only  straight  cracks  can  be 
analyzed.  In  the  present  work,  the  pressure  loading  on  the  local  zone  is  neglected.  As  this  is  a  small  region, 
it  is  good  enough  approximation.  A  better  approach  would  be  to  use  a  different  local  zone  that  can  take  into 
account  bending  also. 

In  the  intermediate  model,  the  adhesive  is  discretized  and  modeled  as  beams.  This  results  in  lack  of 
compatibility  between  the  patch  and  the  sheet.  As  the  discretization  at  this  level  is  not  very  dense  for  the 
adhesive,  the  error  due  to  this  lack  of  compatibility  is  thought  to  be  smedl.  A  possible  remedy  for  this  is  the 
usage  of  an  element  as  the  one  used  in  the  local  model. 

As  the  traction  conditions  are  being  transfered  from  intermediate  to  local,  extra  kinematic  boundary  con¬ 
ditions  have  to  be  applied.  In  the  cases  studied,  the  line  of  the  crack  is  an  axis  of  symmetry  and  therefore  this 
line  could  be  fixed.  But  such  an  approach  might  not  be  applicable  in  other  cases.  The  effect  of  constraining 
other  lines  has  not  been  studied.  A  possible  remedy  is  to  transfer  kinematic  conditions  from  intermediate 
to  local  and  use  a  finite  element  alternating  method  as  in  Wang  and  Atluri  (1995)  which  ca  handle  mixed 
boundary  conditions  on  the  boundary  of  the  local  zone.  Also,  only  in-plane  stresses  are  transfered,  as  a  2-D 
local  analysis  is  used,  the  out-of-plane  stresses  are  neglected.  Again  this  can  be  remedied  if  a  bending  finite 
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element  alternating  method  is  used.  The  transfer  of  kinematic  conditions  from  the  global  to  intermediate 
between  different  models  (beams  -  shells)  could  lead  to  stress  mismatch  at  the  boundary  of  the  models.  The 
error  accumulation  in  such  a  procedure  is  not  known. 

This  analysis  does  not  account  for  the  thermal  stresses  induced  during  manufacture  and  operation.  These 
stresses  will  reduce  the  patch  efficiency.  Fredell  (1994)  studied  the  influence  of  these  stresses  and  has 
demonstrated  that  the  efficiency  reduction  due  to  these  is  severe  and  that  it  cannot  be  ignored.  The  analysis 
done  here  has  to  be  modified  to  take  into  account  the  thermal  effects.  This  analysis  does  not  take  into  account 
the  peel  (normal)  stresses. 
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2.  FORTRAN  code  for  the  alternating  method  for  bending  of  plates 

The  finite  element  part  of  the  code  is  a  modified  portion  of  that  appearing  in  The  Finite  Element  Method  by 
Zienkiewicz  and  Taylor  [7].  Included  with  the  code  are  several  example  problems  that  demonstrate  the 
plane  elasticity  elements  and  the  plate  elements.  The  plate  elements  are  the  DKQ  elements  formulated  in 
“Evaluation  of  a  New  Quadrilateral  Thin  Plate  Bending  Element,”  by  Batoz  and  Tahar  [1].  Explanation  of 
the  examples  is  included  in  Appendix  D.  The  input  mesh  format  is  included  in  this  same  appendix. 

NOTE:  In  the  entire  code,  subroutines  used  from  the  2D  alternating  method  are  in  all  capitals.  Any 
changes  made  to  these  are  in  lower-case  letters.  Therefore,  all  changes  made  by  this  author  are  in  lower¬ 
case  letters,  including  the  initial  finite  element  routines. 

What  is  the  status  of  the  code? 

1.  Using  macro  statements,  plane  elasticity  and  plate  problems  may  be  solved 

2.  The  alternating  method  for  2D  plane  elasticity  was  implemented,  but  never  was  fully  functional.  TJie 
subroutines  were  taken  from  a  2D  code  from  Venkat  Nagaswamy. 

3.  The  alternating  method  for  the  bending  of  thin  plates  is  very  close  to  being  fully  functional.  The  G- 
matrix  which  is  multiplied  by  the  coefficient  matrix  to  get  the  nodal  reactions  seems  to  be  in  order.  For 
the  example,  symmetric  problem  (2x2cracb.inp)  the  G-matrix  appears  to  account  for  the  symmetry  and 
the  boundary  conditions  at  the  end  of  the  first  iteration  appear  to  be  of  the  correct  order  and 
symmetrical.  The  problem  seems  to  be  resolving  the  mesh  using  these  new  boundary  conditions.  For 
some  reason,  the  loads  are  very  small,  but  the  returned  moments  at  the  crack  location  are  too  large.  The 
stress  intensity  factors  then  blow  up  instead  of  converging  upon  iteration.  It  appears,  then,  that  the  only 
problem  is  in  the  latter  part  of  the  iteration  process  where  the  mesh  is  resolved  with  the  new  boundary 
conditions. 

4.  Final  output  needs  to  be  adjusted  to  line  up  correctly  and  to  improve  appearance. 


NOTES: 

1.  The  least  squares  method  for  determining  the  coefficients,  used  in  the  alternating  method  for  bending,  is 
correct. 

2.  The  infinite  integral  solution  functions  (second  half  of  Tmatrix  subroutine  -  located  in  BNDMOM.FOR) 
have  been  verified  using  Mathematica,  The  solution  of  the  entire  infinite  integral  equations  (first  half 
of  Tmatrix  subroutine)  is  not  verified  because  there  is  no  exact  solution,  or  way  of  approximating  the 
exact  solution,  in  some  other  analytic  method. 

3.  Parts  of  the  code  that  deal  with  anisotropic  materials  are  commented.  These  were  to  implemented  later 
when  patching  was  part  of  the  process. 

4.  An  interesting  concept  is  to  implement  the  work  of  Murthy  [3],  which  would  adjust  the  stress  intensity 
factors  to  account  for  crack  closure.  This  is  advantageous  using  the  alternating  method  because  the 
moments  at  the  crack  are  known.  If  s  been  awhile  since  this  paper  was  last  read,  but  the  concept  was 
quite  interesting  at  the  time. 
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Appendix  A 


L  Verification  of  Part  (1)  of  [2] 

A.  The  solution  of  the  transverse  displacement  is  shown  to  be  of  the  form: 

(cj +C2Jc)e 

The  constants  and  Cj  are  solved  for  using  the  given  boundary  conditions.  This  verifies  that 
the  solution  is: 

B.  Given  the  above  form  of  the  displacement  (21),  the  other  seven  functions  are  verified,  which 
are  Equations  (22)-(29). 
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-( - 
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-(x  Z)  -  i  y  Z 
E  fz  1  Z 


2  Dm  (1  -  nu)  (3  +  nu)  P 


-) 


-(— - 

(3  +  nu)  P 

^=-Dm(D[D[D[w,x]  ,x]  ,x] -Z"2  (2-nu)  D[w,x]) 

-  (x  Z)  -  i  y  Z  -  (x  Z)  -  i  y  Z 


3  E 


-  (Dm  (- 


fz  Z  E 
-  + 


fz  Z  (1  +  nu  -  (1  -  nu)  X  Z 


2  Dm  (3  +  nu)  P 


-  (x  Z)  -  i  y  Z 
2  E  fz  E 

(2  -  nu)  Z  ( - - - +  - 

2  Dm  (3  +  nu)  P  Z 
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E  fz  (-1  +  nu)  X  Z 

2  (3  +  nu)  P 

vy=-i  Dm  Z  (Z'^2  w  -  (2-nu)  D[DCwrX],x]) 

-  (x  Z)  -  i  y  Z 


2  Dm  (1  -  nu)  (3  +  nu)  P 
-  (x  Z)  -  i  y  Z 


fz  (1  +  nu  -  (1  -  n 


2  Dm  (1  -  nu)  (3  +  nu)  P  Z 


(E 


fz  (1  +  nu  -  (1  -  nu)  X  Z) ) 


-  (Dm  i  Z  (- 


2  Dm  (1  -  nu)  (3  +  nu)  P 
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G(p)=  ^(l)£”jo«P)sin5d^- Jo«p)sin^ld{. 

If  wc  make  use  of  equation  (2.1.14)  we  find  that  G(p)  =0  when  p  >  1  so 
that  the  equation  (3.2.2)  is  automatically  satisfied  by  a  function  ^({)  of  the 
form  (3.5.1).  Also  wc  note  that 

To  find  the  appropriate  form  of  the  function  we  substitute  from  equa¬ 
tion  (3.5.1)  into  the  integral  defining  the  function 


to  find  that 


F(p)  =  J“r‘^(o^o(?p)d^ 

[  r_m_.  o^p<i 


F(p)  =  J  Jo  V(P  '  ) 

r‘  <^wdi 
Uo  v(p^  -~^y 


p>  1. 


(3.5.5) 


(3.5.6) 


(3.5.7) 


Since  equation  (3.2.1)  is  equivalent  to  the  equation  F(p)=/(p),  0^p<l 
we  see  from  equation  (3.5.6)  that  is  the  solution  of  the  integral  equation 


Jo  V(P^  “ 


=/(p).  0gp<l 


(3.5.8) 


from  which  it  follows  by  the  result  expressed  by  the  equations  (2.3.7)  that  if 
/(p)  is  continuously  differentiable  in  [0, 1] 


^  r 


p/(p)dp 

-p^y 


(3.5.9) 


Substituting  from  equation  (3.5.9)  and  integrating  by  parts,  we  find  that 


at?;  =  |{  cos^  ....... 

which  can  easily  be  transformed  to  give  the  solution 

It  is  easily  verified  that  when  /(p)  =  1,  .0  g  p  <  1.  this  formula  leads  to  the 
expression  (3.1.6)  for  A({). 


u/(u)du 


+  (J  sinitdt 


Jo  v/(‘*  -  “*) 
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In  some  applications  the  quantity 


pG(p)dp 


(3.5.11) 


has  a  physical  significance.  If  we  substitute  from  equation  (3.5.4)  into  equation 
(3.5.11)  we  have 

r'  Pdp  .  f  .  f 

e  =  1 7^)  -  //“'■J.  77^- 

The  first  of  these  two  integrals  can  be  evaluated  immediately;  to  evaluate 
the  second  we  interchange  the  order  of  the  integrations,  carry  out  the  p-in- 
tegration  and  then  perform  the  t-integration  by  using  the  rule  for  integration 
by  parts.  In  this  way  we  find  that 

For  example,  if  /(p)  is  the  polynomial 


(3.5.12) 


fip)  =  I  a,p' 

r-O 


then  it  is  easily  calculated  that 


„  nir  +  i) 

.V'nTTT)- 


(3.5.13) 


(3.5.14) 


We  saw  in  §1.5  how  it  is  possible  to  reduce  the  second  basic  problem  to 
the  first  one  in  a  simple  way,  but  because  it  occurs  so  often  in  applications 
in  physics  and  mechanics  we  shall  show  how  to  derive  (by  an  elementary 
method)  the  solution  of  the  pair  of  dual  integral  equations 


jJmVoiMdi  =  /(p). 

0gp<l 

(3.5.15) 

J^*lA(Wo({p)d{  =  0. 

p>  1. 

(3.5.16) 

In  this  case  we  make  the  representation 

HO  -  x(0sin{/dl, 

z(0)  =  0. 

(3.5.17) 

It  is  obvious  from  equations  (2.1.14)  that  the  form  (3.5.17)  satisfies  the  equa¬ 
tion  (3.5.15)  automatically.  If  we  integrate  by  parts  we  find  that 


it  I  » 
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=  J"jo«p)d{[-x(l)cos{+£x'(Ocos{ld< 

=  J'x'(‘)d/  J*Jo(p{)cos{ld{-x(l)J^  Jo(p«cos{d{. 

Making  use  of  the  result  expressed  by  equations  (2.1.13)  we  see  that  the  integral 
equation  determining  x(0  “ 


r 


x'(»)dt 


=/(p) .  0  ^  p  <  1 . 


Jo  \/(P^  “  **) 

This  is  exactly  the  same  as  the  equation  (3.5.8)  for  <f>(t)  and  so  has  solution 

2  d  f'  p/(p)dp 

If  we  now  integrate  with  respect  to  t  and  use  the  fact  that  x(0)  “0,  we  have 


2  r'  P/(P)dp 

^  J.  J(i‘  -"Fr 


(3.5.18) 


Inserting  this  expression  into  equation  (3.5.17)  we  see  that  the  dual  integral 
equations  (3.5.15)  and  (3.5.16)  have  the  solution 

In  particular,  if  the  free  term  in  equation  (3.5.15)  is  /(p)  •»  1,  0  ^  p  <  1, 

3.6.  Methods  based  on  the  Integral  Representation  of  Harmonic  Functions 

We  saw  in  the  last  section  that  a  solution  of  the  first  basic  problem  is  ob¬ 
tained  if  we  substitute  the  form  (3.5.1)  into  the  harmonic  function  (3.1.2) 
provided  that  is  given  by  equation  (3.5.9).  If  we  make  this  substitution 
then  by  the  analysis  of  the  last  section 

K(p,z)“J  e"«*Jo({p)dij|^  </>(0cos{tdr  (3.6.1) 

is  a  harmonic  function  in  the  half-space  2^0  with  the  property  that 

ei--- 


(3.6.2) 


j/^bic^-1>  <^  -  /\C 


£/®»a\  *<iji  ^•  \  ^ 


/  ^ 

1  ' 

'i(<^  x(A 

a\  -  0 

leT 

Mi,).-  f 

,  1 

\  PA)  COS 
"  / 

U  ^t 

A(t')  -- 

H 

1.': 

PA)  5 

rr 

r 

-iMi), 

-  ft 

OC--  0i'A 

j, .  ( 

r,,  f 

<^“-  TT''^ 


S.-  51 


G(fV 


^^(0  J.tif) 
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Xix)  =  '  A(u)Joixu)du, 


then  we  find  that  dual  integral  equations  of  the  kind 

=  F{x),  Xi(x)  =  G(x) 


(4.4.10) 


9iW  =  t(x),  Xiix)  =  C{x)  (4.4.11) 

arise  in  the  discussion  of  axisymmetric  problems  concerning  the  distribution 
of  stress  in  the  vicinity  of  penny-shaped  cracks  in  infinite  elastic  bodies. 

If  we  substitute  the  values  a  =  v  =  0  into  equation  (4.2.21)  we  find  that 

..X  2  r’  •  j  r'  Tf(T)dT  2.  r*  TG(T)dT 


2u  r‘  ^  f®  TC(t)dT 

"  1  J, 


Similarly  from  equations  (4.3.13)  and  (4.3.14)  we  derive  the  expressions 

uG*(u)du  2x  V(x^  -  1)  f  >  tVd  -  t2)F(t)dr 
^  n  dx  Jj  ^(x^-tt^)  71  Jo  ^  • 


(4.4.13) 


^  2  r  ^  F*(u)du 
“  71  J, 


tG(t)dt 


x^)  71^(1  - X^)  J,  (/2-:c2)^(,2_l)* 


where  F*(u),  C*(u)  are  given  by  the  equations 


F*(u)  =  r  C*(w)= _ ^  f”  5G(5)d5 


(4.4.14) 


(4.4.15) 


In  this  case  it  is  readily  shown  that 


2;:£  xx^(x)dx  =  4  J  xV(l  -  x^)f(x)dx  +  4  J "  4.4.16) 


where 


2f(0  =  2^.(1, 1;  i:  r^),  r>i. 


(4.4.17) 


4.5.  Dual  Integral  Equations  with  Trigonometrical  Kernels 


If  we  write 


/•  00 

(f)(x)  =  J  W~M(w)COSXMdu 
X(.x)  =  J  j4(w)cosx«du 


(4.5.1) 


(4.5.2) 
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then  the  dual  integral  equations 


J  u”‘/4(u)cosxudM  =  F(x),  xe/,, 
J  /4(u)cosxudM  =  0,  X6/2, 


(4.5.3) 


(4.5.4) 


which  arise  in  the  analysis  of  mixed  boundary  value  problems  m  the  plane 
(and  which  are  the  two-dimensional  analogue  of  the  first  pair  of  dual  equa¬ 
tions  considered  in  the  last  section),  can  be  written  in  components  as 


<^,(x)=F(x),  X2(^)=0- 


(4.5.5) 


These  equations  are  of  Titchmarsh  type  with  a  =  /?  =  0,  v  =  -• i.  This  is 

a  particular  case  which  is  not  covered  by  the  Titchmarsh-Busbndge  solution 
(4  22  2)  so  that  special  methods  have  to  be  derived  for  the  soluUon.  Recently 
Srivastava  (1963)  has  considered  the  equations  (4.5.)  with  x(x)  defined  by  equa- 
tion  (4.5.2)  but  with  <f>(x)  defined  by  the  equation 

^(x)  =  J  ‘  ^(w)  cos  XU  du 

but  the  cases  fc=0.  1  are  the  only  ones  which  arc  of  physical  mterest  and  we 
shall  therefore  consider  these  cases  by  the  method  of  Sneddon  (1962c).  ^ 

As  a  special  case  of  the  Wcber-Schaf  heitlin  discontinuous  integral  (2.1.20) 

we  find  that 

J^"j,(u)cosxudu  =  1.  xe/,;  u-'J.(u) cos xudu  =0.  xe/, 

and  that ,  if  ^  is  a  positive  integer, 

J  u"V2,(u)cosxudu  =  (2q)"’ 

J  J2,(u)cosxudu  =0,  xelj 

where  g,  denotes  the  Jacobi  polynomial  defined  by  equation  (2.1.21).  Hence 
we  see  that 

A(u)  =  aouJi(u)  +  21  qa^2,(u) 


wiU  be  a  solution  of  the  pair  of  equations  (4.5.5)  provided  that  the  constants 
are  so  chosen  that 
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Fix)  =  flo  +  f  >■•*).  xelt. 

^•1 


(4.5.7) 


From  equations  (2.6.4)  and  (2.6.5)  we  see  that  the  orthogonality  relation  for 
the  polynomials  5  9  can  be  written  in  the  form 

/•i  ^  ^#0, 

u-ni-«r^5,(0,i,H)5,(0,i,w)du  = 

Jo  ^  =  0, 

where  <5^,  denotes  the  Kroneckcr  delta.  Using  this  relation  and  the  formula 

5,(0, i,«)  =  2F,(-g,9;i;M)=  -m)’"*] 

(Dy  =  d  /du)  we  find  that 


and  that  if  ^  g  1 


=  J_  r*u-*(l  -u)-*F(u*)ciu 
Tt  Jo 


so  that 


(4.5.8) 


2  r‘  F(x)dx 

“  IT  i  v(^’ 

"  rQ)r(g+i)  Jo  9^1-  (4.5.9) 


Since  the  integral 


is  divergent  and  the  integral 


(u) cos  wdu 


J  J2,(«)cosudu  (^^1) 

is  convergent  we  see  from  equation  (4.5.6)  that  if  we  impose  the  additional 
requirement*  on  our  solution  i4(M)  that  x(^)  n^tist  remain  finite  as  x-^  1  — , 

•  In  contact  problems  in  two-dimensional  elasticity  this  condition  is  usually  impos¬ 
ed  by  physical  considerations. 


4.5] 
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then  the  function  F(x)  must  be  such  that  Oo  =0.  wc  see  from  equation 
(4.5.8)  that  this  is  equivalent  to  the  condition 


If  x>  1 


r  F(x)dx  ^ 

Jo 

I  U"‘ Ji,(H)COSXMdu  =  ~2^ 


(4.5.10) 


r  =  [x+y(.x^-l)]- 


so  that  we  obtain  the  formula 


(4.5.11) 


(p2(x)  =  I  (-r)*fl, 
^-1 


(4.5.12) 


where  the  are  given  by  equation  (4.5.9)  and  r  by  equation  (4.5.11).  Similarly, 
by  .making  use  of  the  result 

J  J2^(u)cosxudu  =  (1—  x2)''*5,(0,i,x2),  xeli 


we  find  that 


Xi(x)  -  ^2)  ^ 


Chong  (1953)  considered  the  case  in  which  F(x)  is  a  polynomial 
F(x)  =  f  c„.x". 


(4.5.14) 


If  we  substitute  this  expression  into  equation  (4.5.10)  we  find  that  the  con¬ 
dition  which  must  be  satisfied  if  ^(x)  is  to  remain  finite  as  .x  -*  1  —  is 


£  Hi"  tMc  =  0. 

.To  r(i«  +  1)  " 


(4.5.15) 


This  is  the  criterion  derived  otherwise  by  Chong. 

If  F(x)  is  a  polynomial  containing  only  even  powers  of  x,  the  results  are 
much  simpler.  Suppose,  for  example,  that 


F(X)=  I  C2„X^ 

n-0 


(4.5.16) 


then  it  is  easily  shown  that 
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[  V  (qgn), 

d;f(«*)=  j  "•»  ”• 

0,  (9>n) 

and  hence  from  equation  (4.5.9)  that  a,  =  0  if  ^  >  n  and  that 

.fcg'x  (,&”>■  (‘•.S.I7) 

’  V’T  »-o  n!r(23  +  n  +  l) 

It  should  also  be  noted  that  if  XiW  =  0  and  ^,(x)  is  a  polynomial  of  degree 
m  in  X*  then 

y  (x)  =  (4.5.18) 

where  n„(x’)  denotes  a  polynomial  of  degree  m  in  x 2,  and  ^2(*)  is  a  poly¬ 
nomial  of  degree  2m  in  9  =  [x  +  “  i)]  • 

If  F(x)  =cosnKx  then  it  is  easily  shown  that 


(4.5.17) 


(4.5.18) 


and  it  follows  from  equation  (4.5.9)  that 


2  r*  cosnTTDdi^  ,  .  . 

;  J. 


so  that  the  solution  corresponding  to  F(x)  *=  cos  httx  is 


Hence  if 


<^4(u)  =  «Ji(u)Jo(n7r)  +  4  2  qJ2^(n7t)J2fl(w) . 


F(x)  —  T.  b„  cos  MJTX.  X  e  / 1 

ir-O 


(4.5.19) 


the  solution  of  the  dual  integral  equations  (4.5.5)  is 


A(u)  =  uJi(u)  I  b„Joinn)  +  411  (4.5.20) 

fr«0  n-O?-! 

If  x(x)  is  to  remain  finite  as  x  —  1  - ,  the  coefficients  b,  must  satisfy  the 
linear  relation 
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flo  +  2  aJoO"^)  =  0- 


(4.5.21) 


The  solution  (4.5.20)  was  first  derived  by  Fredricks  (1958). 

In  a  similar  way  it  can  be  shown  that  if  F(0)  =  0,  the  solution  of  the  dual 

integral  equations 


J  i4(«)sin.xudu  =  F(.x),  xeli 

^  00 

J  y4(u)sinxud«  =0,  xelj 


can  be  expressed  in  the  form 


I  (2q  +  l)V2,+  i(«)» 

f  «o 


(4.5.22) 


(4.5.23) 


J  u***(l  -  u)’"*D;[u'*F(«*)]du  (4.5.2 


“*■  r(i)r(9  +  i)  Jo" 

(cf.  Sneddon,  1963c,  §5).  From  this  it  can  readily  be  deduced  that 

eo 

F(x)  =  1  6,sinn>rx,  xef, 

»-  I 

then  the  solution  of  the  pair  of  equations  (4.5.22)  is 

A(u)  =  £  £  (2q  +  l)  M2,+  •(""V 2,+ 1(“) 

II-I 

in  agreement  with  Fredricks’  result  (Fredricks,  1958). 

If  wc  write 


(4.5.25) 


(4.5.26) 


if,(.x)  =  J^"u 


i/l(u)cos  XU  du 


(4.5.27) 


and  define  x(x)  by  equation  (4.5.2)  then  the  equations  (4.5.5)  are  of  the  type 
we  encounter  in  discussing  crack  problems  in  the  two-dimensional  theory 
of  elasticity.  An  elementary  solution  of  these  equations  can  be  found  by  a 
method  similar  to  that  used  in  §3.5.  If  we  represent  A{u)  in  terms  of  a  function 
g(r)  by  the  formula 


A(u)  =  g(0.^o(»“)d< 


(4.5.28) 
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then  substituting  in  equation  (4.5.2)  and  interchanging  the  order  of  the  in¬ 
tegrations  we  find  that 


X(x)  =  f  g(0df  f  Jo(tu)cosxudu 
Jo  Jo 


and  using  the  result  (2.1.13)  we  see  that,  in  this  case, 


VO' 


XrW  =  0. 


(4.5.29) 


Similarly  if  we  write 


=  -  f' 

dx  Jo 


A(u)  sin  XU  du , 


substitute  the  form  (4.5.28)  for  A(u),  interchange  the  order  of  the  integrations 
and  make  use  of  equation  (2.1.14)  we  find  that 


-si 


ir(0dt 

V(^'-'')’ 


fO)dt 


lo 


(4.5.30) 


It  follows  immediately  that  (4.5.28)  will  give  a  solution  of  the  dual  integral 
equations  (4.5.5)  provided  g(0  is  a  solution  of  the  integral  equation 


-  r 

dx  Jo 


=  F(x),  xel, 


It  is  easily  shown  by  means  of  equations  (2.3.7)  that 


_2r  r 

'  Jo 


F(w)dw 

^(l2  _  ^2)- 


(4.5.31) 


The  solution  of  the  dual  integral  equations  is  given  by  equations  (4.5.28) 
and  (4.5.31)  and  the  auxiliary  functions  XiW*  by  equations  (4.5.29), 
(4.5.30)  and  (4.5.31).  It  should  also  be  noted  that  the  first  equation  of  the  pair 
(4.5.29)  leads  to  the  result 


XiWd* 


and  that  substitution  of  the  fonn  (4.5.31)  for  g(t)  into  this  equation  gives 
the  formula 


f  =  r  V(1  - 

Jo  Jo 
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A  solution  of  this  type  has  been  derived  by  Tranter  (1963a)  for  the  pair  of 
equations  (4.5.3)  and  (4.5.4).  If  wc  integrate  equation  (4.5.4)  with  respect 
to  X  we  obtain  the  equation 


J  u^^i4(u)sinxudu  =  C,  x€/2, 

where  C  is  a  constant.  Now  it  is  well-known  that 

lim  r  ii”^.r4(M)sinxud«  =  i7r.4(+0) 

x-^co  Jo 

and  if  the  integral  <f>{x)  defined  by  equation  (4.5.1)  is  to  exist  wc  must  have 
y4(+  0)  =  0  so  that  C  =  0  and  equation  (4.5.4)  is  equivalent  to 

J  u"“^i4(u)sinxudM  =  0.  (4.5.32) 

If  now  we  multiply  equation  (4.5.3)  by  (p^  —  x*)"^  and  integrate  with  respect 
to  X  from  0  to  p  wc  find,  on  making  use  of  the  relation 

cosuxdx  1  T  /  \ 

obtained  by  applying  Hankel’s  inversion  theorem  to  equations  (2.1.13),  that 

0S(>< 

Similarly,  if  wc  multiply  equation  (4.5.32)  by  (x^  —  p^)’*  and  integrate  with 
respect  to  x  from  p  to  oo ,  we  find  that 

j  u*" * A(u)Jo(pu)du  =  0 ,  p  >  1 . 

Applying  the  Hankcl  inversion  to  the  last  two  equations  wc  find  that  the  func¬ 
tion  A{u)  is  given  by  the  equation 

(4.5.33) 

Similarly  to  solve  the  pair  of  equations  (4.5.22)  we  diflferentiate  the  first 
equation  with  respect  to  x  to  obtain  the  equivalent  relation 

j*  i4(u)cosxiidw  =  xe/i. 

Proceeding  in  the  same  way  as  in  the  previous  case  we  can  then  show  that 
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A{u)J  q{pu)6u  =  H(\ 


-’/.'sS 


(x)dx 


x^) 


from  which  it  follows  by  an  application  of  Hankel’s  inversion  theorem  that 

AW  =  I  ,4.5.34) 


4.6.  Dual  Integral  Equations  with  Hankel  Kernel  and  Arbitrary  Weight  Function 

In  this  section  we  shall  consider  the  solution  of  dual  integral  equations 
of  the  type 


Xiix)  =  G(x)  (4.6.1) 

where  the  functions  F{x)  and  C{x)  are  prescribed  in  the  intervals  and 
respectively  and  the  functions  4>ix)  and  x(x)  are  defined  in  terms  of  ^he  known 
function  k{u)  and  the  unknown  function  A(u)  by  means  of  the  equations 


(pix)  =  u~^*[l  +  k(«)]/l(u)7^(xu)d«. 


A{u)J  ^{xu)du. 


(4.6.2) 


If  we  write  A{u)  =  u^iu)  and  define  two  functions  f{x)  and  g{x)  by  the 
equations 


fix)  =  S^,_.,2,{[1  +  k(i/)]<A(H);x},  g{x)  =  S^,^o'P(x)  (4.6.3) 

then  the  equations  (4.6.1)  are  easily  shown  to  be  equivalent  to  the  pair 

/,(x)  =  2^“;c“^'f(;c),  g2(x)  =  Gix).  (4.6.4) 


4.6.1.  Reduction  to  a  Fredholm  Equation 
We  can  reduce  the  problem  of  solving  the  dual  integral  equations  (4.6.4) 
to  that  of  solving  a  Fredholm  integral  equation  of  the  second  kind  by  using 
a  method  similar  to  that  employed  in  §4.2.4  above.  Putting  /?  =  0  in  equation 
(4.2.45)  we  make  the  representation 


'I'M  =  5+,._./i(u)  (4.6.5) 

of  the  unknown  function  {l/(u)  in  terms  of  a  new  unknown  function  h(u). 
Substituting  from  equation  (4.6.5)  into  equations  (4.6.3)  we  find  that 


and  that 


fix)  —  S.jv-a.2»^iv,-«/>(^)  +  Siy-t.2a{Kx)S^y  ^ghix)} 
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foi 


S  oc  a  <3  w 


A  ^  ^  ^ 


r  i’>^ 


A  i.  -  ■  \  r  c  ;  <-v  ^ 


A(uV-  -o(A^t 


9C(.) 


C  '  '  }a 


r^s  'Xu. 


o  ■' 


r  >  ■^ 


'A- 


1/  y  y 


X,  ^  c 


0(i)  ~  Jx  •^'*'^  ■^u  c<t 


•'  i" '  rl  -  \  '■ 
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Appendix  C 


I.  When  Equations  (47)  and  (49)  are  substituted  into  the  even  and  odd  parts  of  Equations  (30)-(37)  and 
(38)-(45),  respectively,  the  resulting  equations  are  of  the  form  shown  in  these  handwritten  notes.  You 
will  realize  that  both  parts  of  each  function  are  not  fully  written  out  because  the  minor  differences  in  the 
equations  may  be  accounted  for  by  observation. 

II.  The  infinite  integrals  are  solved  analytically  by  a  method  similar  to  Sneddon,  pp.  496-497  [5].  The 
derivations  and  symbols  are  slightly  different  than  those  used  in  [2]. 

The  following  is  equivalent  to  Equation  (75),  (142),  and  (143)  in  [5]: 


C(a,A:,y)=  cos%y)J di. 

(75) 

S;(a,A:,y)=  sin(^)/„(a^^'^"cf^ 

(142) 

C  /■S'”(<3t,x,y)  =  Z;  (a,x  +  iy) 

(143) 

NOTE:  The  Mathematica  output  is  utilized  in  the  handwritten  portion.  The  notation  is  the  same. 
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Z=l/(a"2+w^2)^(l/2) 

1 


2  2 
Sqrt [a  +  w  ] 

Integrate [Z,wl 

2  2 

Log  [w  +  Sqrt  [a  +  w  ]  ] 

ZOlsw  +  Sqrt  [a'^2+w^2] 

2  2 
w  +  Sqrt [a  +  w  ] 

w=x+I  y 
X  +  I  y 
ZOl 

2  2 

X  +  Sqrt [a  +  (x  +  I  y)  ]  +  I  y 
ComplexExpand [Expand [ZOl] ] 

22  2  2  22  1/4 

x+{4x  y  +(a  +  x  -  y)) 

2  2  2 
Arg[a  +x  +  2Ixy-y] 

Cos  [ - - - - — ■]  + 

2 

22  2  2  221/4 

I(y+(4x  y  +(a  +  x  -  y)) 

2  2  2 
Arg[a  +  x  +  2Ixy-y] 

Sin[ - ]) 

2 

ro={ (r  Cos [O] +R^ (1/2) Cos [Fe/2] ) ^2+ (r  Sin[0]+R^ (1/2) Sin [Fe/2] ) ^2) ^ (1/2) 
Fe  2 

Sqrt  [  (Sqrt  [R]  Cos  [ — •]  +  r  Cos  [O]  )  + 

2 

Fe  2 

(Sqrt[R]  Sin[ — •]  +  r  Sin[0])  ] 

2 

Simplify [ro] 

2  Fe 

Sqrt  [r  +  R  +  2  r  Sqrt  [R]  Cos  [ - O]  ] 

2 

Zo2=(a^2+(x+l  y)^2)^(-l/2) 
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2  2 
Sqrt  [a  +  (x  +  I  y)  ] 
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AltsoW 


Comp lexBxp and [Expand [Zo2] ] 
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(4x  y  +(a  +  x  -  y)) 
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Arg[a  +  x  +  2Ixy-y] 
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2 


22  2  2  22  1/4 

(4x  y  +{a  +  x  -  y)) 

Z03=(x+I  y)/(a'‘2+(x+l  y)*2)"(3/2) 

X  +  I  y 


2  2  3/2 

(a  +  (x  +  I  y)  ) 

ComplexExpand [Expand [Z03]  ] 

2  2  2 
3  Arg[a  +  x  +  2Ixy-y] 

X  Cos  [ - - - 

2 

22  2  2  22  3/4 

(4x  y  +{a  +  x  -  y)) 


2  2  2 
3  Arg[a  +  x  +  2Ixy-y] 
y  Sin  [ - — — •] 
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22  2  2  22  3/4 

{4x  y  +(a  +  x  -  y)) 

2  2  2 
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(4x  y  +(a  +  x  -  y)) 
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3  Arg[a  +  x  +  2Ixy-y] 

X  Sin[ - - •] 


22  2  2  223/4 

(4x  y  +(a  +  x  -  y)) 

Z04=-l/(a"2+(x+I  y)"2)"{3/2)+3(x+I  y) "2/ (a"2+ (x+I  y)"2)"(5/2) 

2 

2  2  -  (3/2)  3  (x  +  I  y) 

-  (a  +  (x  +  I  y)  )  + - 

2  2  5/2 

(a  +  (x  +  1  y)  ) 
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ComplexExpand [Together [ExpandAll tZ04] ]  ] 

2  2  2 
2  2  2  5  Arg[a  +  x  +  2Ixy-y] 

(-a  +  2  X  -  2  y  )  Cos  [- - - - ■] 


22  2  2  22  5/4 

(4x  y  +(a  +  x  -  y)) 


2 

5  Arg  [a 
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2 

+  X  +2 
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I  X  y 


22  2  2  22  5/4 

(4x  y  +(a  +  x  -  Y  )  ) 


+ 


2  2  2 
5  Arg [a  +x  +2Ixy-y] 

4  X  y  Cos  [ - - - ■] 


I  (- 


22  2  2  22  5/4 

{4x  y  +{a  +x  -  y  )  ) 


2  2  2 
22  2  5  Arg [a  +x  +2Ixy-y] 

(-a  +  2  X  -  2  y  )  Sin[- — — - •] 


-) 


22  2  2  22  5/4 

(4x  y  +{a  +  x  -  y  )  ) 

Z05=-3/2  G  (x+I  y)/(a^2+{x+I  y) ^2) * (5/2) +5/2  G  (x+I  y) "3/ (a^2+ (x+I  y) ^2) 


3 

-3  G  (x  +  I  y)  5  G  (x  +  I  y) 


2  2  5/2  2  2  7/2 

2  (a  +  (x  +  I  y)  )  2  (a  +  (x  +  I  y)  ) 

Together [Z05] 

2  3  2  2  2 

(-3  a  GX  +  2GX  -  3Ia  Gy+6IGx  y-6Gxy  - 

3  22  27/2 

2IGy)/(2(a  +  x  +  2Ixy-y)  ) 


3 
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ComplexExpand [Together [Z05]  ] 

2  3  2 
((-3a  GX+2GX  -6Gxy) 

2  2  2 

7  Arg[a  +  x  +  2Ixy-y] 
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2 

22  2  2  22  7/4 

(2  (4  X  y  +  (a  +  X  -  y  )  )  )  + 

2  2  3 

((-3a  Gy+6Gx  y-2Gy) 

2  2  2 

7  Arg[a  +  x  +  2Ixy-y] 

Sin  [ — - — - ■]  )  / 

2 

22  2  2  22  7/4 

(2  (4  X  y  +  (a  +  X  -  y  )  )  )  + 

2  2  3 

I  (((-3a  Gy+6Gx  y-2Gy) 
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2 
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2 


22  2  2  22  7/4 

(2  (4  X  y  +  (a  +  X  -  y  )  )  )  ) 

Z06=  (p^2+a^2)  ^  (-5/2)  *G5*LegendreP  [4,p/  (p^2+a''2)  ^  (1/2)  ] 

4  2  2  4 

G5  (3  a  -  24  a  p  +  8  p  ) 


2  2  9/2 

8  (a  +  p  ) 

p=x+I  y 
X  +  I  y 
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2 
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15  a  X  Sin[ - ■] 

2 

- — )  - 

22  2  2  22  9/4 

(4x  y  +(a  +x  -  y)) 

2  2  2 
9  Arg [a  +x  +2Ixy-y] 

-15  a  X  Cos  [ - ] 

2 

8  X  y  ( - - - — - 

22  2  2  22  9/4 

{4x  y  +(a  +  x  -  y)) 

2  2  2 

9  Arg [a  +x  +2Ixy-y] 

15  a  y  Sin[ - ■] 

2 

- - - )) 

22  2  2  22  9/4 

(4x  y  +(a  +  x  -  y)) 
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Appendix  D 


I.  A  copy  of  the  user  instructions  from  Zienkiewicz  [7]  is  included. 

II.  Additional  macro  commands  are  also  given. 


III.  Example  Problems 
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Additional  Macro  Statements 


TABLE  1 

CRACK  DEFINITION 


CRACK  DEFINITION  -  Format(il0,4fl0.4)  -  must  immediately  follow  a  CRAC  macro. 

The  crack  type  is  as  follows:  kedge(ncrac)->  crack  type  (l=edge;  0  or  2=embedded;  3  S3mim) 


Example:  2x2crac.inp 

Columns 

Descriptions 

Variable 

1-10 

Crack  Type 

kedge(ncrac) 

10-20 

X-Coord.  - 1 

xyck(l,l,ncrac) 

20-30 

Y-Coord.  -  1 

xyck(2,l,ncrac) 

30-40 

X-Coord.  -  2 

xyck(l,2,ncrac) 

40-50 

Y-Coord.  -  2 

xyck(2,2,ncrac) 

TABLE  2 
SYMMETRY 

SYMMETRY  -  Forraat(4iI0)  -  must  immediately  follow  a  SYMM  macro. 

The  variables  SYMIA  and  SYMIB  give  any  two  nodes  on  the  X  axis  of  symmetry.  Likewise  the 
variables  SYM2A  and  SYM2B  give  those  on  the  Y  axis  of  symmetry.  Note  that  the  axes  of  symmetry  can 
only  be  parallel  to  X  and/or  Y  axes  but  need  not  be  the  X  and  Y  axis. 


Columns 

Descriptions 

Variable 

1-10 

SYMIA 

10-20 

SYM2A 

20-30 

SYMIB 

30-40 

SYM2B 

TABLE  3 

DISTRIBUTED  LOAD 

DISTRIBUTED  LOAD 

-  Fonnat(3il0,4fl0.0),  Format(5fl0.0)  -  must  immediately  follow  a  DIST  macro. 

Two  lines  that  give  the  element  N,  the  generator  increment  NG,  the  Gaussian  Rule  LS,  and  the  remaining 

parts  give  the  loads  at  each  gauss  point.  If  only  one  load  is  given,  it  is  assumed 

to  be  evenly  distributed 

over  the  element.  Note: 

TWO  LINES  MUST  ALWAYS  BE  USED  Example: 

PLATESS.INP 

Columns 

Descriptions 

Variable 

1-10 

Element  Number 

N 

10-20 

Generator  Increment 

NG 

20-30 

Gaussian  Rule 

LS 

30-40 

Gauss  pt#l 

ld(l) 

40-50 

Gauss  pt  #2 

ld(2) 

50-60 

Gauss  pt  #3 

ld(3) 

60-70 

Gauss  pt  #4 

ld(4) 

10-20 

Gauss  pt  #5 

ld(5) 

20-30 

Gauss  pt  #6 

ld(6) 

30-40 

Gauss  pt  #7 

ld(7) 

40-50 

Gauss  pt  #8 

ld(8) 

60-70 

Gauss  pt  #9 

ld(9) 

TABLE  4 

BOUNDARY  MOMENT 


BOUNDARY  MOMENT  -  Fonnat(4iI0,9fl0.0),  Format(5fl0.0)  -  must  immediately  follow  a  MOME 
macro.  Distributed  moment  along  boundary  edge  at  specified  gauss  points.  Example:  2x2cracb.inp 

Columns  Descriptions  Variable 


1-10 

Node  Number  1 

N1 

10-20 

Node  Number  2 

N2 

20-30 

Generator  Increment 

NG 

30-40 

Gaussian  Rule 

LS 

40-50 

Mx  -  Pt.  #1 

Id(l) 

50-60 

Mx  -  Pt.  #2 

ld(2) 

60-70 

Mx  -  Pt.  #3 

Id(3) 

10-20 

My-Pt.#l 

ld(4) 

20-30 

My  -  Pt.  #2 

ld(5) 

30-40 

My  -  Pt.  #3 

ld(6) 

40-50 

Mxy  -  Pt.  #1 

ld(7) 

60-70 

Mxy  -  Pt,  #2 

ld(8) 

70-80 

Mxy  -  Pt.  #3 

ld(9) 

TABLE  6 

POINT  MOMENT  LOAD 

POINT  MOMENT  LOAD  - 

Format(2fl0.5,fl0.2)  -  must  immediately 

follow  a  MOMP  macro. 

Columns 

Descriptions 

Variable 

1-10 

X  Coord. 

XC 

10-20 

Y  Coord, 

YC 

20-30 

Mx 

LD(1) 

30-40 

My 

LD(2) 

TABLE? 

DISTRIBUTED  SHEAR  LOAD 

DISTRIBUTED  SHEAR  LOAD  -  Fonnat(4il0,6fl0.0)  -  must  immediately  follow  a  DISQ  macro. 
Distributed  shear  load  between  nodes  N1  and  N2  at  up  to  6  gauss  points.  Example:  plateldq.inp 


Columns 

Descriptions 

Variable 

1-10 

Node  Number  1 

N1 

10-20 

Node  Number  2 

N2 

20-30 

Generator  Increment 

NG 

30-40 

Gaussian  Rule 

LS 

40-50 

Gauss  pt#l 

id(i) 

50-60 

Gauss  pt  #2 

ld(2) 

60-70 

Gauss  pt  #3 

ld(3) 

70-80 

Gauss  pt  #4 

ld(4) 

80-90 

Gauss  pt  #5 

ld(5) 

90-100 

Gauss  pt  #6 

ld(6) 

TABLE  8 

POINT  SHEAR  LOAD 

POINT  SHEAR  LOAD-  Format(2fl0.5,fl0.2)-  must  immediately  follow  a  POIN  macro. 
Example:  plateld2.inp 

Columns  Descriptions  Variable 


1-10 

10-20 

20-30 


X  Coord. 

Y  Coord. 
Shear  Load 


XC 

YC 

LD 


Example  Problems 


All  the  example  problems  that  follow  are  done  on  the  same  simple  mesh  shown  in  Figure  1.  All  element 
lengths  are  1 .0,  therefore,  each  side  has  a  length  of  2.0 

Ya 


7  8  9 


1  2  3 


Figure  1.  Simple  mesh  used  in  all  examples  (all  elements  lengths  =  1.0) 


Description 

Force  in  X  and  Y  direction  at  Node  9 
My  moments  between  1-3  and  7-9 
2D  alternating  method  example 
Alternating  method  for  bending  example 
Simply  supported  plate  with  distributed  shear  load 
Cantilever  example  with  boundary  shear  load 
Point  shear  load  example 


File  Name 

examp  le2.inp 

2x2plate.inp 

2x2crac.inp 

2x2cracb.mp 

platess.inp 

plateldq.inp 

plateld2.inp 


example2.inp 


feap**  Example  Problem 
9, 4, 1,2, 2,4,0 


coord 
1.1, 0,0 
3,0,2,0 
4,1,0,1 
6,0,2,! 
7,1,0,2 
9,0,2,2 

elem 

1.1.1.2.5.4.1 

3.1.4.5.8.7.1 


boun 
1.1.1, -1 
3,0,0, 1 
4,0,1,0 
7,0,1,0 


fore 

9,0,5,5 

mate 

1,1 

100,.3,1.0,3,3,1 

end 

macr 

tang 

form 

solv 

disp 

stre 

end 

stop 


2x2plate.inp 


feap**  Example  Problem 
9,4, 1^,3, 4,0 


coord 
1,1, 0,0 
3,0,2,0 
4,1,0,1 
6,0,2,1 
7,1, 0,2 
9,0, 2,2 

elem 

1.1.1.2.5.4.1 

3.1.4.5.8.7.1 


boun 
1,0,1,0,0 
2,0,1,0,0 
3,0,1,0,0 
7,0,1,0,0 
8,0, 1,0,0 
9,0, 1,0,0 


mome 

1,2,0,3,0,0,0 

10,10,10,0,0,0 

2,3,0,3,0,0,0 

10,10,10,0,0,0 

9,8,0,3,0,0,0 

10,10,10,0,0,0 

8,7,0,3,0,0,0 

10,10,10,0,0,0 


crac 

2,.8,1.,1.2,1. 

mate 

U 

1000000,.3,0.04,2,2,1 

end 

macr 

tang 

form 

solv 

disp 

stre,node 

end 

stop 


2x2crac.mp 


feap**  Example  Problem 
9,4, 1,2, 2, 4,0 


coord 
1,1, 0,0 
3,0,2,0 
4, 1,0,1 
6,0,2,! 
7,1, 0,2 
9,0,2,2 


elem 

1.1.1.2.5.4.1 

3.1.4.5.8.7.1 

boun 
1,1, 1,-1 
3, 0,0,1 


fore 
7,0,0,5 
8,0,0,10 
9, 0,0, 5 

crac 

2,.8,1.,1.2,1. 

mate 

1,1 

100, .3, 0.04, 2,2,1 


end 

mao- 

tang 

form 

solv 

disp 

stre,node 

afem,memb 

end 


2x2cracb.inp 


feap**  Example  Problem 
9,4,1,2,3A0 


coord 
1,1, 0,0 
3, 0,2,0 
4, 1,0,1 
6,0,2, 1 
7,1, 0,2 
9, 0,2, 2 


elem 

1.1.1.2.5.4.1 

3.1.4.5.8.7.1 


boun 
1,0,1,0,0 
2,0, 1,0,0 
3,0, 1,0,0 
7, 0,1, 0,0 
8,0, 1,0,0 
9,0, 1,0,0 

mome 
1, 2,0,3, 0,0,0 
10,10,10,0,0,0 
2,3,0,3,0,0,0 
10,10,10,0,0,0 
9,8,0, 3, 0,0,0 
10,10,10,0,0,0 
8,7,0,3,0,0,0 
10,10,10,0,0,0 


crac 

2,.8,1.,1.2,1. 

mate 

1,2 

1000000,.3, 0.04,2,2,1 

end 

macr 

tang 

form 

solv 

disp 

stre,node 

afem,bend 

end 


platess.inp 


feap**  Example  Problem 
9.4, 1^,3, 4,0 


coord 
1,1,0,0 
3,0A0 
4,1.0,1 
6,0,2,! 
7,1, 0,2 
9,0,2,2 


elem 

1,1,1A5,4,1 

3,1,4,5,8,7.1 


boon 
1,0,1,0,0 
2,0,1,0,0 
3,0,1,0,0 
4,0,1,0,0 
6,0, 1,0,0 
7,0, 1,0,0 
8,0, 1,0,0 
9.0,1,0,0 

dist 

1,0,3,100 

2,0,3,100 

3,0,3,100 

4,0,3,100 

mate 

1,2 

10000000,-3,0.04,3,2,1 


end 

macr 

tang 

form 

solv 

disp 

stre,node 

end 

stop 


piateldq.inp 


feap**  Example  Problem 
9, 4, 1,2, 3,4,0 


coord 
1.1,0, 0 
3, 0,2,0 
4,1,0,1 
6, 0,2,1 
7,1,0,2 
9,0,2,2 


elem 

1.1.1.2.5.4.1 

3.1.4. 5. 8. 7.1 


boun 
1.0,1,1,1 
4,0, 1,1,1 
7,0,1, 1,1 


disq 

3,6,0,3,100. 

9,6,0,3,100. 

mate 

la 

10000000,-3, 0.04,2,2,1 


end 

macr 

tang 

form 

solv 

disp 

stre,node 

end 

stop 


plateld2.inp 


feap**  Example  Problem 
9,4, 1,2,3, 4,0 


coord 
1,1, 0,0 
3,0,2,0 
4,1,0,1 
6,0^1 
7,1, 0,2 
9,0,2,2 


elem 

1.1.1.2.5.4.1 

3.1.4.5.8.7.1 

boun 

1,0,1,1,1 

2,0,1,1,1 

3,0,1,1,1 

4,0,1,1,1 

7,0,1,1,1 


pom 

2.,2.,100. 

mate 

10000000,. 3, 0.04^A1 

end 

macr 

tang 

form 

solv 

disp 

stre^iode 

end 

stop 


Appendix  C 


I.  When  Equations  (47)  and  (49)  are  substituted  into  the  even  and  odd  parts  of  Equations  (30)-(37)  and 
(38)-(45),  respectively,  the  resulting  equations  are  of  the  form  shown  in  these  handwritten  notes.  You 
will  realize  that  both  parts  of  each  function  are.  not  fully  written  out  because  the  minor  differences  in  the 
equations  may  be  accounted  for  by  observation. 

II.  The  infinite  integrals  are  solved  analytically  by  a  method  similar  to  Sneddon,  pp.  496-497  [5].  The 
derivations  and  symbols  are  slightly  different  than  those  used  in  [2]. 

The  following  is  equivalent  to  Equation  (75),  (142),  and  (143)  in  [5]: 


C,7(a,x,>')= 

(75) 

(142) 

C  (a,x,y)- i  S„"  (a,x,y)  =  z;  (a.x  +  iy) 

(143) 

NOTE:  The  Mathematica  output  is  utilized  in  the  handwritten  portion.  The  notation  is  the  same. 


AltsoI4 


Z=l/(a"2+w^2)"(l/2) 

1 


2  2 
Sqrt [a  +  w  ] 

Integrate [Z,w] 

2  2 

Log [w  +  Sqrt [a  +  w  ]] 

201=w  +  Sqrt [a^2+w^2] 

2  2 
w  +  Sqrt [a  +  w  ] 

wssx+I  y 
X  +  I  y 
ZOl 

2  2 

X  +  Sqrt [a  +  (x  +  I  y)  ]  +  I  y 
Con^lexExpand  [Expand  [ZOl]  ] 

22  2  2  22  1/4 

x+(4x  y  +{a  +  x  -  y)) 

2  2  2  ^ 
Arg[a  +  x  +  2Ixy-y] 

Cos  [ - — - ]  + 

2 

22  2  2  221/4 

I(y+(4x  y  +(a  +  x  -  y)) 

2  2  2 
Arg[a  +  x+  2Ixy-y] 

Sin[ - - - ]  ) 

2 

ro= ( (r  Cos [0] +R^ (1/2) Cos [Fe/2] ) ^2+  (r  Sin [O]  +R^ (1/2) Sin [Fe/2] ) "2) ^  (1/2) 
Fe  2 

Sqrt [ (Sqrt [R]  Cos [ — ]  +  r  Cos [0] )  + 

2 

Fe  2 

(Sqrt[R]  Sin[ — ]  +  r  Sin[0])  ] 

2 

Simplify [ro] 

2  Fe 

Sqrt [r  +  R  +  2  r  Sqrt [R]  Cos  [ - O] ] 

2 

Zo2=(a"2+(x+I  y)^2)^(-l/2) 

1 


2  2 
Sqrt [a  +  (x  +  I  y)  ] 


1 


AltsoM 


CoiKpIexExp2uid  [Expsuad  [Zo2]  ] 

2  2  2 
Arg[a  +  x  +  2Ixy-y] 
Cos  E - - 


22  2  2  221/4 

(4x  y  +(a  +x  -  y)) 


2  2  2 
Arg [a  +x  +2Ixy-y] 

I  sin[ - - - •] 


22  2  2  22  1/4 

(4x  y  +(a  +  x  -  y)) 

Z03=(x+I  y)/(a^2+(x+I  y)^2)^(3/2) 

X  +  I  y 

2  2  3/2 

(a  +  (x  +  I  y)  ) 

Comp lexExp and  [Expand  [Z 03]  ] 

2  2  2 
3  Arg [a  +x  +2Ixy-y] 

X  Cos  [— - — - ] 

2 

22  2  2  22  3/4 

(4x  y  +(a  +x  -  y)) 

2  2  2 
3  Arg [a  +x  +2Ixy-y] 

y  Sin[ - — - ■] 

2 

22  2  2  22  3/4 

(4x  y  +{a  +  x  -  y)) 


2  2  2 

3  Arg [a  +x  +2Ixy-y] 

y  Cos  [ - 

2 

I  ( — - - — ; - 

22  2  2  223/4 

(4x  y  +(a  +  x  -  y)) 

2  2  2 

3  Arg [a  +x  +2Ixy-y] 

X  Sin[ - — - ■] 


22  2  2  223/4 

(4x  y  +{a  +  x  -  y)) 

Z04=-l/(a^2+(x+I  y)^2)^(3/2)+3(x+I  y) ^2/ (a*2+ (x+I  y)^2)‘(5/2) 

2 

2  2  -  (3/2)  3  (x  +  I  y) 

-  (a  +  (x  +  I  y)  )  +  — - - - - 

2  2  5/2 

(a  +  (x  +  I  y)  ) 


2 


Altsol4 


Comp lexExp and [Together [ExpandAll [Z041 ] ] 

2  2  2 
2225  Arg[a  +x+2Ixy-y] 

(-a  +  2x  -  2y)  Cos  [ - — - 3 

2 


22  2  2  22  5/4 

(4xy+(a+x-y)) 

2  2  2 
5  Arg [a+x+2Ixy-yl 

4  X  y  Sin  [ - - - ] 

2 

- - -  + 

22  2  2  22  5/4 

(4x  y  +(a  +  x  -  y)) 

2  2  2 
5  Arg [a+x+2Ixy-y] 

4  X  y  Cos  [ - - — - - 3 

2 

I  ( - — — -  - 

22  2  2  22  5/4 

(4x  y  +(a  +  x  -  y  )  ) 

2  2  2 

2225  Arg [a  +x  +2Ixy-y] 

(-a  +  2x  -  2y)  SinC - — - 3 

2 

- - - ) 

22  2  2  22  5/4 

(4  X  y  +  (a  +  X  -  y  )  ) 

Z05=-3/2  G  (x+I  y)/(a^2+(x+I  y) ^2) * (5/2) +5/2  G  (x+I  y) ^3/ (a^2+ (x+I  y)*2)' 


-3  G  (x  +  I  y) 


5  G  (x  +  I  y) 


2  2  5/2  2  2  7/2 

2  (a  +  (x  +  I  y)  )  2  (a  +  (x  +  I  y)  ) 

Together [Z05] 

2  3  2  2  2 

(-3a  GX+2GX  -3Ia  Gy+6IGx  y-6Gxy 

3  22  27/2 

2IGy)/(2(a  +  x  +  2Ixy-y)  ) 


3 


Altsoi4 


Con^ilexExpamd  [Together  [Z05]  ] 

2  3  2 
((-3a  GX+2GX  -GGxy) 

2  2  2 

7  Arg[a  +  x  +  2Ixy-yl 

Cos  [ - - - ■] )  / 

2 

22  2  2  22  7/4 

(2  (4  X  y  +  (a  +  X  -  y  )  )  )  + 

2  2  3 

((-3  a  Gy+6Gx  y-2Gy) 

2  2  2 

7  Arg[a  +  x  +  2Ixy-y] 

Sin  [— - - - ■] )  / 

2 

22  2  2  227/4 

(2  (4  X  y  +  (a  +  X  -  y  )  )  )  + 

2  2  3 

I  (((-3  a  Gy+6Gx  y-2Gy) 

2  2  2 

7  Arg[a  +  x  +  2Ixy-y] 

Cos  [— - — - 3  )  / 

2 

22  2  2  227/4 

(2  (4  X  y  +  (a  +  X  -  y  )  )  )  - 

2  3  2 

((-3  a  Gx  +  2Gx  -  6Gxy) 

2  2  2 

7  Arg[a  +x  +2Ixy-y] 

Sin  [— — - 3 )  / 

2 

22  2  2  22  7/4 

(2  (4  X  y  +  (a  +  X  -  y  )  )  ) ) 

Z06=  (p'^Z+a^S)  (-5/2)  *G5*LegendreP  [4,p/  (p''2+a''2)  (1/2)  ] 

4  2  2  4 

G5  (3  a  -  24  a  p  +  8  p  ) 

2  2  9/2 

8  (a  +  p  ) 

p=x+I  y 
X  +  I  y 


Altsol4 


CoznplexExpand  [Z06] 

4  2  2  2 

(G5  (3  a  -  24  a  (x  -  y  )  + 

4  2  2  4 

8(x  -  6x  y  +  y)) 


Cos  [~ 


9  Arg[a  +  (x  +  I 


2 

y)  ] 


-])  / 


2  2  9/2 

(8  Abs  [a  +  (x  +  I  y)  ]  )  + 

2  3  3 

(G5  (-48  a  xy+8(4x  y-4xy)) 

2  2 
9  Arg  [a  +  (x  +  I  y)  ] 

Sin  [ - ]  )  / 

2 

2  2  9/2 

(8  Abs  [a  +  (x  +  I  y)  ]  )  + 

2  3  3 

I  ( (G5  (-48  a  xy+8(4x  y-4xy)) 

2  2 

9  Arg  [a  +  (x  +  I  y)  ] 

Cos  [ - — — ] )  / 

2 

2  2  9/2 

(8  Abs  [a  +  (x  +  I  y)  ]  )  - 

4  2  2  2 
(G5  (3  a  -  24  a  (x  -  y  )  + 

4  2  2  4 

8(x  -  6x  y  +y)) 

2  2 

9  Arg  [a  +  (x  +  I  y)  ] 

Sin  E - - - 3  )  / 

2 

2  2  9/2 

(8  Abs  [a  +  (x  +  I  y)  ]  ) ) 


5 


AltsoMb 


Zll=  (  (a"2+ (x+I  y)  "2)  ^.(1/2)  -  (x+I  y)  ) /a 
2  2 

-X  +  Sqrt [a  +  (x  +  I  y)  ]  -  I  y 


Comp lexExp and  [Expand  [Zll]  ] 

2  2  2 
2  2  2  2  2  2  1/4  Arg  [a  +  x  +  2Ixy-y] 

(4x  y  +{a  +  x  -  y))  Cos  [ - ^ - ] 

X  2 

-(-)  + - - - ^ ^ -  + 

a  ^ 

2  2  2 
2  2  2  2  2  2  1/4  Arg  [a  +  x  +  2Ixy-y] 

(4x  y  +(a  +  x  -  y))  Sin[ - : 

y  2 

I  (-(-)  + - — - - - 

a  3. 

Z12=l/a- (x+I  y)  / (a* (a^2+ (x+I  y)^2)^(l/2)) 

1  X  +  I  y 

a  2  2 

a  Sqrt [a  +  (x  +  I  y)  ] 

ComplexExpand  [Expand  [Z12  ]  3 

2  2  2 
Arg [a  +x  +2Ixy-y] 

X  Cos  [ - - - - — '] 

1  2 _ 

a  2  2  2  2  2  2  l/4 

a  (4  X  y  +  (a  +  x  -  y  )  ) 

2  2  2 
Arg [a  +x  +2Ixy-y] 

y  sin[ - ■] 

2 

- — -  + 

22  2  2  22  1/4 

a(4x  y  +(a  +  x  -  y)) 

2  2  2 
Arg  [a  +  x  +  2Ixy-y] 

y  Cos  [ - - - ■] 

2 

I  (-( - - - ^  + 

22  2  2  22  1/4 

a(4x  y  +{a  +  x  -  y)) 

2  2  2 
Arg  [a  +  x  +  2Ixy-y] 

X  Sin[ - - - ] 

2 

- - - ) 

22  2  2  22  1/4 

a  (4  X  y  +  (a  +  x  -  y  )  ) 

Z13=a/(a"2+(x+I  y)"2)"(3/2) 


2  2  3/2 

(a  +  (x  +  I  y)  ) 


Altsol4b 


Comp lexExp and  [Expand  [Z13  ]  3 
2  2 

3  Arg[a  +  x  +  2  I  x  y 

a  Cos  [ - - - - - 

2 


2 


y 


] 


-3 


2  2  2  2 
(4  X  y  +  (a  +  X 

2  2 

3  Arg  [a  +  x 
I  a  Sin[ - 


+ 


2 


2  2  3/4 

y  )  ) 

2 

2  I  X  y  -  y  ] 
- - — ^1 


22  2  2  223/4 

(4x  y  +(a  +  x  -  y)) 

Z14=3 (x+I  y) ( (x+I  y)+(a^2+(x+I  y) ^2) ^ (1/2 ) ) / (a  2+(x+I  y) 
a/((x+I  y)+(a"2+(x+I  y)^2)"(l/2)) 

3  a  (x  +  I  y) 


2  2  5/2 

(a  +  (x  +  I  y)  ) 

ComplexExpand  [Expsind  [Z14]  ] 

2  2  2 
5  Arg [a  +x  +2Ixy-y] 

3  a  X  Cos  [ - - - 

2 

22  2  2  22  5/4 

{4x  y  +(a  +  x  -  y)) 

2  2  2 
5  Arg [a  +x  +2Ixy-y] 

3  a  y  Sin[ - - - - 

2 

22  2  2  225/4 

(4x  y  +(a  +  x  -  y)) 


3  a  y  Cos  [- 


2  2 
5  Arg  [a  +  x  + 


I  (- 


2  2  2  2 
(4  X  y  +  (a  +  X 


2 

2  I  X  y  -  y  ] 
- ] 


2  2  5/4 

y  )  ) 


2  2  2  . 
5  Arg  [a  +  x  +  2Ixy*-y]^ 
3  a  X  Sin[ - — - - - 3 


-) 


22  2  2  225/4 

(4x  y  +{a  +  X  -  y)) 


2)"(5/2)* 


1 


Altsol4c 


P15=IjegendreP  [3  ,  -1,  z] 

2  2 
- ( (1  -  5  z  )  Sqrt [1  -  z  ] ) 

8 

z=(x+I  y)/(a"2+(x+I  y)"2)"(l/2) 

X  +  I  y 

2  2 
Sgrt [a  +  (x  +  I  y)  ] 

Z15=24/(a^2+(x+I  y)^2)^2*Pl5 

2  2 
5  (X  +  I  y)  (x  +  I  y) 

-3  (1 - -)  Sqrt  [1 - - - •] 


2  2 
a  +  (x  +  I  y) 


2  2 
a  +  (x  +  I  y) 


2  2  2 
(a  +  (x  +  I  y)  ) 

PowerExpand [%] 

2  2  2 
3a  (-a  +4x  +8Ixy-4y) 

22  27/2 

(a  +  x  +  2Ixy-y) 

ComplexExpand [Sin^lify [%] ] 

2  2  2 
2227  Arg[a  +x  +2Ixy-y] 

3  a  (-a  +  4  X  -  4  y  )  Cos  [ - — - ’] 


22  2  2  227/4 

(4x  y  +{a  +  x  -  y  )  ) 

2  2  2 
7  Arg[a  +  x  +  2Ixy-y] 
24  a  X  y  Sin[ - - - 3 


22  2  2  227/4 

(4x  y  +(a  +  x  -  y)) 

2  2  2 
7  Arg[a  +  x  +  2Ixy-y] 
24  a  X  y  Cos  [- - - - 3 


I  (- 


22  2  2  22  7/4 

(4x  y  +(a  +  x  -  y)) 

2  2  2 
2227  Arg[a  +x  +2Ixy-y] 
3a  (-a  +  4x  -  4y)  Sin[ - - - 3 


-) 


22  22  22  7/4 

(4x  y  +(a  +  x  -  y)) 


1 


AltsoMc 


AitsoWc 


ComplexExpand [%] 

2  2  2 
9  Arg[a  +  x  +  2Ixy-y] 

-15  a  y  Cos  [ - — - — - 3 

2 

8  X  y  { - - — - -  + 

22  2  2  22  9/4 

{4xy+(a+x-y)) 

2  2  2 
9  Arg[a  +  x  +  2Ixy-y] 

15  a  X  SinC— — — — - 3 

2 

- — - ^ - )  + 

22  2  2  22  9/4 

(4x  y  +(a  +  x  -  y)) 

2  2  2 
(3  a  -  4x  +  4y)  ( 

2  2  2 
9  Arg[a  +  x  +  2Ixy-y] 

-15  a  X  Cos  [- - — - 3 

2 

22  2  2  22  9/4 

{4x  y  +(a  +  x  -  y)) 

2  2  2 
9  Arg[a  +  x  +  2Ixy-y] 

15  a  y  Sin[ - - - ■] 

2 

- - - — - —)  ■  + 

22  2  2  22  9/4 

(4x  y  +(a  +  x  -  y)) 

2  2  2 
I  ((3  a  -  4x  +  4y) 

2  2  2 

9  Arg[a  *fx  +  2Ixy-y] 

-15  a  y  Cos  E - - - 3 

2 

( - - -  + 

22  2  2  22  9/4 

(4x  y  +(a  +  x  -  y)) 

2  2  2 

9  ArgCa  +  x  +  2Ixy-y] 

15  a  X  Sin[ — — - ■] 

2 

- - - - 

22  2  2  22  9/4 

(4x  y  +{a  +  x  -  y)) 

2  2  2 
9  Arg[a  +x  +  2Ixy-y] 

-15  a  X  Cos  [ - ^ - 3 

2 

8  X  y  (  . . — - - 

22  2  2  229/4 

(4x  y  +{a  +x  -  y)) 

2  2  2 
9  Arg[a  +  x  +  2Ixy-y] 

15  a  y  Sin[ - * - “3 

2 


22  2  2  22  9/4 

(4x  y  +(a  +  x  -  y)) 


3 


Appendix  D 


I.  A  copy  of  the  user  instructions  from  Zienkiewicz  [7]  is  included. 

II.  Additional  macro  commands  are  also  given. 


III.  Example  Problems 


Appendix  E 


Copy  of  “A  finite  element  alternating  approach  fro  the  bending  analysis  of  thin  cracked  plates”  by  Chen 

[2] 


J^ffiuuionalJiwnuthtf  tnuturi’ Ut:  ,  , 
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n  .  .1  .  i  d  ni  ne  llieorv  Ihe  nnalviienl  snliiiimi  for  an  iiitinie  iliin  plalc  cnniainini:  a  eraek 

ni accurtie  nmtL  ciLtUL  ,r  i  .m.i  iin-  innuoncc  of  Ihc  ucomcinc  Inmiulitrics  on  the 

u^r  «..innK'  tracks  The  imcrjctjoii  elTeci  amoni:  eiaeks  amt  the  iniMitntt  oi  im. 

ill"  .  . . . . . . .  . . .  “ 

iajonstniltf  llte  valiJiiy  of  approach. 


I.  Introduction 

The  analysis  of  a  plaic  coniaioing  cracks  subjccicil  lo  bc.uling  lypc  loavlings  has  been 
svesiigaied  exlcnsively  during  ihe  past  iwo  decades  because  of  ns  sipniliaince  in  various 
tneinelring  applications.  One  lo  ibe  existence  of  cracks,  the  serviee  strength  o  a  cracked 
pla'tewill  be  diminished  because  of  the  singular  cbaraciertslies  of  the  stress  held 
•rack  tin  Fll  Hcncc,  the  accurate  ticicrminalion  of  bcinling  stress  intensity  factors,  w  nc  t 
'jdicale  the  severity  of  cracks,  shouhl  receive  tmieh  aiieniion  for  the  safe  design  of  plate 

“unliirthe  ease  for  the  plate  snbjeeted  to  ,n-plate  loadings,  the  prediction  ..f  bending  stress 
latensity  factors  is  dinereiit.  In  tlie  literature,  the  analytical  solnitons  lor  bending  cracked  plates 
rre  available  only  for  limiied  problems  will,  simple  geomeiry  and  boundary  conditions. 
Warns  [I]  obtained  the  stress  distribution  in  ihe  vicinity  ol  the  crack  tip  for  plates  loaded  in 
bending,  and  found  that  the  stress  was  prop.>riional  to  the  inverse  of  the  .siiuare  root  of  the 
adial  distance  from  the  crack  tip.  Sib.e.  al,  1 2 1  n.sed  Williams-  results  to  dehne  symme.ne  and 
,„,i.syinmciric  stress  intensity  factors  for  cracked  plate  problems  involving  bending  h.ads. 
Incorporated  with  Reissner’s  transverse  shear  theory,  the  elleet  of  plate  Ibiekness  on  the  analysis 
o( cracked  plates  in  bending  was  then  siudicil  by  several  researchers  |  3  71.  To  deal  with  moie 
realistic  problems,  however,  numerical  approaches  such  as  the  linite  element  metluuls  appear  to 
^imperative.  Wilson  and  Thompson  IS)  iirst  calculated  the  beiuling  stress  intensity  laetor  f..r 
j  rectangular  plate  containing  a  center  crack,  subiccted  to  pure  cylindrical  bending  by  a 
conventional  linite  element  method.  To  aceonnl  for  the  singularity  ol  the  stress  held  near  the 
crack  lip,  various  types  of  singular  cicincnis  ;uul  liiulc  clcmcni  iiunlcU  were  then  tlevdopul  .uu) 
employed  [9  15].  Sosa  and  r.isehen  fl(i|  proposed  an  iiuliicct  wa\  to  compute  the  bending 

stress  intensity  factor  via  a  paih-indepeiuleni  iniegial.  ,  .  , 

Among  lhe.se  surveys,  the  accuracy  ol  the  lieiuling  stress  iniensiiy  laciors  ohiamed  largely 
depends  on  mesh  si/e  and  expensive  comin.ici  time  is  consiimal  especially  loi  solving  the 
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prohicm  with  multiple  cracks  l  o  compensate  (i>r  these  sluuicomin^s.  an  cnident  tcdiniquc  that 
comhincs  tltc  atlvantaues  ol  (inite  element  metlKul  and  Schwar/s  altemating  method  [17]  is 
thus  sliulicd  in  this  work.  As  is  known,  the  so*calleil  tiniie  element  aliernalinj!  method  has  been 
successfully  applied  in  solvinu  ihiee  dimensnmal  iracture  problems  (ISj  and  two  dimensional 
fracture  problems  \\^)  21  |.  I'or  the  analysis  of  cracked  plates  i.n  bendiiiii.  however,  no  work  has 
been  done.  The  aim  of  this  work  is  thus  to  cMciul  the  linilc  clement  alternating  method  to  deal 
with  the  beiuling  of  a  thin  plate  with  smele  tu  multiple  cracks.  liascd  on  the  classical  plate 
theory  [22]  aiul  the  solution  oluamed  l\v  Sncildon  ami  I  liiott  \  2}\  for  the  plate  with  a  crack 
under  symmetric  normal  pressure,  the  analytical  stdution  lor  an  infinite  thin  plate  containing  a 
crack  for  which  the  crack  sui laces  are  subjected  to  arbitrary  symmetric  bending  loadings 
represented  by  a  pidynonual  o\  appiopriate  order  is  (irst  derived.  An  iterative  superposition 
procedure  is  then  used  to  satisfy  traction  and  displacement  boundary  conditions  on  the  external 
boundaries  of  the  plate  repeatedly  until  the  variation  of  bending  stress  intensity  factors  on  each 
crack  is  negligible. 

I'o  demonstrate  the  .applicability  and  versatilitv  of  the  new  technii|ue.  several  numerical 
c.\amples  concerning  thin  plates  with  single  i>r  multiple  ciacks  subicctcd  to  various  types  of 
loatling  and  bouiularv  ei^mlitions  aie  analyml  Hie  mllucnce  of  the  boundaries  on  the 
computation  of  bending  siiess  intensity  factois  aiul  the  interaction  cnect  among  cracks  arc  also 
studied  in  detail 


2.  Analytical  soluliun  for  a  eraekeil  thin  filatc  in  tieiuling 


based  on  the  classical  plate  theor>  \Z2\,  in  the  absence  of  distributed  lateral  loading, 
governing  equation  for  a  thin  plate  can  be  written  as 

"  .VVVV  '  Ik 

Tiic  relations  among  the  lateral  dcllcciion  »v  and  related  variables  arc 

=  -Ofu'.,,  +  nv,J. 
niyy  -  -  />(!  -  v)w 
,1^  =  -  l){\\\yy,  I  U-.Urk 

1/v  ~  ^7(u' I  u 

r\  =  </.  I  I  (2  v)\\  ,,J 

and 

/*,  ^  ‘/v  -1-  niy 


the 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 
(7) 


/>|U  ,,,  I  (2 


(8) 


4 
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. .  ,i„  hcnding  ami  iwisling  huhikiUs.  irans- 

ihcrc  uiiii  Icnpll'  «f  •'"=  " 

Z,T  . . . .  '■  * . .  ''■'“'■'  . . . . 

rti...,. m...  .^o  '*->  “>■ 

nnsform  technique  [24]  as; 


O 

li 

_ 

1 

d.v^ ' 

.1 

m,  = 

-/)(«■•..%»  -  's''''. 

= 

-  ;)(>•>;•„„  -  s-’u). 

»I.X,  = 

;  j(i  -  v)/);"'.,,. 

10  (H)  become 


(9) 

(10) 
(in 
(12) 
(i.i) 

(14) 

(15) 

(K>) 

(17) 


=  i/)^(«\,>  -  s""  '- 

-(2- v);^v.,) 

inJ 

=  iD4((2  -  »•)«•,..  -  s""') 

The  solmion  of  (9)  is  thus  ohtaiitcii  as  follows: 

uiv.4)  =  (<-,  - 

,  .r  ,iu-  virnblc  c  and  can  be  determined  by 

.« . I  . . 

the  given  boundary  eoiulitions.  Ihe  tran.  |.,,nrier  transform  for  these  variables,  one 

LLw  ..o.„ ,  .0, ...  ,.7..  M... 

obtains  the  analytical  solution  eoiisisteni  I 

. . . .  . . 

.1.1  Iliiii  nl  ites  in  bending,  the  solution  to 
To  obtain  Ow arbitrary  edge  bending  load  must  be 
,he  problem  of  a  semi-inlinite  plate  ‘  ,,„,a  „elow. 

ohtaLd.  .see  I'.g.  1-  The  boundary  eoiulitions  foi  tins  pr. 

(i)HU=/(.v)  at  v  =  0.  ,,,  .1,  II  IS  salislied  auioinatieally.)  .1 

(From  the  ci|uilibrium  condition.  )  , 


% 
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(ii)  i  v  =  0 
and 

(iii)  -♦  0 


/  t,t  I  A  siMtii-injirnU-  plalc  stiluctlctl  It*  piifc  Iviuliitj: 
al  V  -  0 


:  jAfa-  '^r/0  , 


By  condition  (lii),  the  constants  <  ,  and  (4  in  (17)  must  vanish,  l  luis.  (17)  can  he  simplified  to: 

Mt\\s=)  =  (c,  +  (  .A)c  (18) 

Carrying  out  llic  roiiricr  transform  on  both  sides  of  condition  (i)  and  (ii).  we  obtain 


m,  =  Hi) 


(19) 


and 


r,  =  0  PO) 

where  iii)  =  Jl  ,  /(,)•)<-•'’’  il.f  .Subsiiiuiini;  (IS)  inio  Kiins.  (10)  ;iikI  (15)  and  comparing  with  tlie 
variable."!  m,  and  r,  obtained  by  (I9|  and  (20).  the  consliinls  c,  and  t  j  are  determined  and  the 
transformed  «(.v.  s)  in  (IS)  can  be  thus  foiinil  as 


>ttv.{)  = 


-  I 


(.1  +  i  )(i  -  v)i>\i 


!  /(s')l(l  +'•)-(!  -  v)U=|vJc-l’'' 


(21) 


Now.  substituting  (21)  into  Ikiiis.  (10)  (16)  and  talcing  the  inverse  Fourier  transform 
respectively,  the  analytical  solution  of  the  bending  semi-inlinite  thin  plate  is  derived  as: 


-  I 


./(<) 


Id  +  >)-(l  -  v)k'|.v]c-'’l''-‘''d4. 


/k)L(.*  i  '■)-  (I  +v)|el,v|c  I''-'  '-'de. 


(22) 

(23) 


L 
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—  { I  —  v).v  . ,  .  , . 

I\(.v.  r)  =  — -r - c  '  c(is(..tliL- 

n(.1 V)  |„ 


and 


rV(A.  y)  = 


7r(3  +  V) 


'•»  II 


v)j\|c  ''sinUr)d;; 


(b)  for  odd  funclion  f^iv) 


utv.  y)  = 


n(3  -f  v)(l  -  y)I) 


‘  /d-l, 


I  V]  (1  »L  v)c  '' sinic y) dJ. 


y)  = 


ff(3  4-  V)  J 


/*:(c)|(3  I  V)  il  I  r)j\|c  '' sin(.:rld:. 


r)  =  1  I  w\|i:  ^'sinijrldc, 

7r(3  f  V)  „ 


y)  = 


I 


q  jx.  y)  ’ 


n(3  +  V)  ^ 

-2  r 


?r(3  +  V)  j 


/:(v)L-2  4(1“  'Iwyv.lc  *'cos{;yldi. 


</r(A,  y)  = 


r.,(.v,  y)  = 


7r(3  4*  V)  J 
-(1  - 


.1  ^  *  '  Ct>s(kry)  tlw. 


n(3 


—  v).\  f  '  •  .  - 

-  /,(«:);“  c  "'smu:y)d: 

4-  v)  J,, 


and 


rv(v,  y) 


_ 1 _ * 

;r(3  +  V)  ^  „ 


/:(';)sl(''  ■  '■)  -*  (1  -  vljA  )c“’'cos|cy)dc. 


In  ihc  above,  llic  fnnclions  /i(^)  aiul  /:(»:)  arc  dclincti  as 


lAi) 


/i(  y)ci>s(Jy)  dy 
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and 


0  «  i, . . .  "  ™  ‘  "•  . . . 

on  crack  sia  faces 

1"J'm  ars,I.)lc.cJ . . 

crack  closure  on  the  compressuiu  side  is  observtil.) 


(i)  m.  =  /(  V)  =  Z  ' 


at  v  -  0. 1  rl  ^  "• 


(jal!ll'  =  0  at  .V  =  0.1  rl  >  i(. 
(•X 

(ill)  I'.  =  0  at  A  =  0 


and 


(iv)  IV  -♦  0  as  A  -  ±  J. 

Due  to  the  symmetry  of  the  beiuliui:  momem  m..  only  the  region  v  »  0  needs  to  be  analy/ed. 

.!■<  "«  to..,,,*  . . .a...  ....  ■>  '■>  »  “ 

solutions  for  even  ii  aiubor  odd  ii.  lespeetivcly 


1(H)  U CVii'M  ct  til. 

(a)  For  even  orilcr  « 

Let  /i(r)=  n  '  <  ..O  npi»M  siihsiitiitine  llic  hoiuulary  o»iulilions  (i)  aiul  |ii)  into  (31) 

aiul  (30).  a  pair  of  dual  iiilcgral  eiiiialions  lesull: 


/,  (C)  cos(ir)clJ  ^ 


ii: 

/ ,  (cjcosicy)  . 


/,(  D 


0 


Define  the  follouin^  paiaincters; 


/j 

a 


(46) 


y 


ViOi- 

./i(v')=  2/.''-/-,(/>) 


aiui 


r  *1*  * 

cos(/ii;)  =  I  II 


where  Jr  ilciioics  ihe  Hcsscl  function  of  llic  (irsi  kimi  wlili  onlci  r  {lure  r  =  -  1).  l•t|u:llion  (46) 
can  be  rcwrillen  as; 


f'' 


The  solution  of  can  be  obtained  Ironi  Sneddon  124)  as: 


2,A'^ 


Juiff) 


r’'’(!  r')‘  '(/it  v)dr 


'f />  f  »'''(!  i/‘)‘  *1  I  </i(  ra)r' ■./,(/M’)ilr  iln  • 


/••.(P)  = 
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Hence,  from  the  acf.niiu.n  as  siaicl  ab..vc.  one  olnains 

■  1  f 


r(  ~  +  2 


■(•  ac  1  r'  ■  1  dr  . 

1. 

■J..  1 

found. 

(b)  For  odd  order  ii 

..  - 1.'. .....  .....i  . . . . . .  . . . 

m«.  (.TOI  >,..l  (Ml. ..  p..i. ..( 

\ 

« J..  V 

I  /j(4)sin(ir)  V  =  0  y  >  >1  J 

Define  ihc  fnneiions  /..(;)  anil  sinl,-'/)  as  follows: 

(Jii'l)  =  .A* 

and 


=  0  ,1'  >  >1 


sin(,»;)  =  ["y-]' 

Equation  (4K)  can  be  rewiilien  as; 

I-'/-'* 

/■■'l/'Ui  ''  '■  ' 

J" 

From  Sneddon  124].  /  .(/O  is  deiiveil  as 

/s  .M  iv  r> 


r’-’H  t-l' -./..ivldr 

\  n  /  I 

r  *  *  *  * 

I,.  ic'-ll  Ii-l' •  >;.(  vKiv' ■./.■(/>itdi- 

Jll  .  “ 


(49) 


i  ••  / 

/ 
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and  ./j(i;)  is  ihus  observcil  as: 


/:(C) 


(I'n 


s 

V 


•  1 

-  (■„ 
) 

./ ,  (wJ)  H  i/J 

r" ' "./  >|f/w  r)<l  V 
n 

/ 


SuhstiUiling  (49)  ituo  l-Apis.  (.VS)  (45),  (lie  itnniyiical  si»lulinns  ft*!  u-.  /fi,,  r,  and 

r,.  due  1()  nu'incnl  Io;ulini!  uilh  tnld  otilcr  n  :uc  thus  obuiiiicd 


3.  Computniion  of  stress  intensity  faeinr 

In  ucncml.  dcfortnalions  near  a  crack  up  occur  in  three  dilVereni  modes,  namely,  opening 
mode  (mode  I).  sliding  mode  (nuule  M)  and  tearing  mode  (mode  111).  The  severity  of  the 
crack-tip  slrcss/displacemenl  licltl  is  ilcsctihed  hy  three  ci>rrespoiHling  parameter.s 
Kyi  which  are  known  as  stress  intensity  facii>rs  ( 25  |.  I'or  a  cracked  plate  in  bending. ‘how¬ 
ever,  the  fracture  behavior  is  describeii  by  iljc  symmetric  and  atui-symmelric  stress  intensity 
factors  as  deiined  by  .Sih  ei  al.  (31.  Since  only  the  symmetric  type  problem  is  considered  in 
this  work,  the  symmetric  bending  sitess  intensity  factor  which  is  actually  the  nuHle  I  stress 
intensity  factor  K,  is  therefore  stiulied.  I  ioni  the  ilelinition  by  ITdogan  [25]  and  the 
relation  between  the  stress  component  a,,  and  moment  ni,  in  plate  theory  [22],  K|  can  be 
expressed  as 


K,  =  lim  ~j^/2(r  -  </)//!, (O.  r). 

.r-»«  ft 

Substituting  the  aforementioned  analytical  solution  of  into  the  above  expression,  and 
performing  tedious  manipulations.  Kt  is  obtained  as: 

(a)  for  even  order  it 

(b)  for  odd  order  n 


k-  -  V  r  r 


I'he  complete  expression  for  K ,  is  thus’ 


'  '.-.r:,,..  Ir  2"|(»/2)!|-  '  2” "  !|  (;i/2)/2 1!!  ’ 
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Substituting  (47)  and  (49)  into  Iit|ns.  (30)  (37)  and  (38)  (45)  respectively,  the  analytical 
solutions  for  u’,  moment  loading  with  order  n  arc 

obtained.  Tlicy  can  be  expressed  in  terms  t)f  the  coellicienls  C'„(ii  =  0,  1.2,...)  in  matrix  form 
as; 


(50) 


{K)  =  [S]{Cl 


(51) 


where  the  column  vectors  \Q\,  \  K  ;  ami  \('\  are  ilelineil  as: 


IQ}  —  f  U’, /lly, /Mj., /ilxv*  */r«  ^  i‘ 1 


and 

(C}  =  [Co,C, . 

([7]^  [5])  arc  the  matrix  of  relation  function  between  the  analytical  solution  vectors  |Q)  and 
(X)  for  u', /Mjr, /Hy, /»,r» * *» *V  ^  '  cocllicient  vector  {Cj.  Both  [T^and  [s]  arc 

functions  of  crack  length.  In  addition,  [T]  is  also  a  function  of  space.  As  a  result,  once  the 
coefficient  vector  {C)  is  determined,  the  vector  {Q]  at  any  location  of  the  infinite  plate  and  the 
bending  stress  intensity  factor  {K  ]  can  be  calculated  from  (.50)  and  (51). 


4.  Finite  element  ulternating  procedure 

The  finite  clement  alternating  procedure  employed  by  ('hen  and  Chang  [19  21]  for  two 
dimensional  problems  is  extcmlcd  io  solve  the  present  cracked  plate  problems  in  bending. 
Accordingly,^  finite  cracked  plate  problems  can  be  superpi^sed  by  finite  domain  problems 
without  cratks  and  infinite  domain  proldems  with  a  crack  subjected  to  suitable  beiuling 
moments  on  crack  surfaces,  l•‘urlhermore.  multiple  crack  problems  can  Ik  .su|Krposed  by  single 
cracked  problems  succes.sively  in  the  similar  manner.  I  he  fast  <lecaying  nature  and  self- 
equilibrating  condition  of  the  residual  bending  moments  on  the  crack  surfaces  and  resultant 
residual  iraction/displacemenl  tm  the  external  boundary  show  the  great  advantage  of  the 
method  in  obtaining  covergenl  solutions,  l  or  completeness,  ihe.se  procedures  arc  described 
briefly  as  follows: 

<l)  Input  geometry  and  boundary  conditions  of  the  bouiuled  crackeil  plate  to  be  analyzed. 

^2)  fivaluatc  the  residual  bending  moment  nt{i  at  the  location  of  fictitious  crack 
k{k  =  1,2, ...» Nr)  by  a  conventional  finite  element  method  ('fhe  crack  need  not  lie  along  clement 
boundaries.).  Nc  is  the  total  number  of  cracks  existing  in  the  plate. 

<3>  Since  the  crack  surfaces  are  free  of  tractions,  reverse  the  sign  of  the  residual  bending 
moment  m{i  and  fit  it  with  a  polynomial  of  appropriate  order.  The  filled  residual  bending 
moment  can  be  expressed  as  where  is  the  polynomial  cocllicieni 

vector  for  crack  k  and  1/,J  is  the  vector  composed  of  (  —  0.  1,  2 . N)  in  which  u*  is  the 


/ 
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half  length  of  the  Ath  craek.  l  o.  bcitc.  aecuraey.  '.c;!'  can  he  ilctcrnnncd  hy  a  least-squares 
incthoil  through  the  ininiini/alion  ol  the  integral. 


(/«(;-  »i"i I'll"  tk  -i.:. 


{cir  =  [i/]:’!«;r. 


l«li= 

J  -wi. 


<4>  Since  the  syinnietiie  hciuling  type  problems  are  solved,  the  inieraetion  cHeet  among 
eneks  is  induced  from  the  bending  moment  m,  only.  I  rom  ICI"  in  step  (M  and  (50).  the  reversed 
residual  inleraetir.n  eoellicient  veeiors  indueal  by  craek  k  at  lictitious  eraek  f(/ =  1.  2.....Wc, 

I  ^  k)  {di„  Kin  be  evaluated.  If  /  >  A.  arid  '..I'u  to  IC’!'  and  compute  the  updated  ’D,.  Thus 
the  cocmeicnl  vector  {C\l  of  the  residual  bending  moment  for  craek  A  needed  to  lie  released 
for  the  next  computation  is  fmuul  to  be  ICi;  =  li'A. .  .  Id'.r,.  Kepeat  all  measures  as  mention^ 
above  until  the  cumulative  residual  bending  moments  on  each  crack  arc  negligible.  As  a  result, 
the  resultant  eoenieient  vector  for  crack  k\CU  =  !f!l’  +  \nl  +  •  f‘>r  the  current  iteration  is 

<5>  Substituting  the  resultant  eoeHieient  vector  for  eraek  AjC  j»  into  (.50)  and  (51).  the 

parameters  iv./«,.».,..../,../,.r,.iv  and  K,  can  be  computed  and  uiidaied. 

<6>  Once  the  variation  of  the  beiuling  stress  intensity  factor  K|  for  all  eraek  lips  is  smaller 
than  a  small  value  say. 

I K, (current  iteration)  <  .. 

K, (current  iteration) 

the  iteration  is  (inished.  Here,  y  is  taken  as  ().(H)5. 

<7>  To  satisfy  loading  conditions  on  ti)e  external  houndary,  the  lesultant  residual  equivalent 
nodal  force  !C?i:  on  the  external  boundary  at  the  alb  element  is  eoinpuled  from  (.50)  by  the  sum 
of  the  contribution  of  each  crack  iyj,A  as 


!(>!?=  I  !y!.r. 


where 
(C}ai  = 
and 


/ 

/ 
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[G],k  dcnoles  Ihc  relation  matrix  and  ,s,  is  the  honiulary  where  loading  conditions  applied.  [A/] 
iipi^ms  the  sluipe  fiinelion.  [l/^i  is  the  irimsformaiinn  matrix  ami  f  7')i  is  the  matrix  of  relation 
(unction  for  the  kiU  crack  as  ilelineii  in  (50).  To  satisfy  the  ilisplacement  Imittttlary  wmlition  ot>  the 
external  bonmlary.  the  resultant  resiilual  notlal  ilispkicemcnts  jiij*  in  the  external  ilisplacement 
boundary  at  node  /( arc  computeil  from  (50)  by  the  sum  of  the  contribution  of  each  crack  !i(|,u  as; 

("}?  =  Z 

1 


where 

{«},»  =  mi.uf'ii 


[ZV  =  [A/],v[r],. 

ts  the  iransformalioii  matrix  aiul  {  ly  is  tldincil  in  (M)). 

(8)  Reverse  the  resultant  residual  ei|uivalcni  luulal  lorcc  lt?It  and  restdlanl  residu.d  luulal 
displacement  \a];  calculated  in  step  (7)  and  cimsidet  those  as  new  li)ading  aiul  displacement 
conditions  acting  on  the  external  boundarie.s.  Repeat  alt  the  iteratit)n  steps  <2>  to  <7>  until  the 
variation  of  bending  stress  intensity  facti)i  K,  on  each  crack  tip  is  negligible. 


5.  Results  and  discussions 

To  verify  the  validity  of  the  new  propcKsetl  method,  several  plate  heiuling  problems  with  single  or 
double  cracks  under  difll-rent  moment  and  displacement  houiulary  condmons  are  solved.  The 
effect  of  interaction  between  cracks  and  the  inlluence  4>r  boundaries  on  the  compulation  of  beiuling 
stress  intensity  factors  are  also  studied.  As  for  finite  element  motlcling,  the  simple  four-node  plate 
element  leach  node  has  four  tiegrees  of  freedom  u  .  u\„.,  n  .,.  and  u  is  employed  in  all  analyses. 

Example  I:  A  stpuire  center  cracked  thin  plate  suhiccied  to  itni/intn  hendintf 

As  seen  in  Fig.  3  and  Fig.  4,  a  square  center  cracked  plate  subjected  to  uniform  bending  moment 
along  the  boundary  at  y  =  ±  //  is  first  solved.  The  geometry  and  its  (inite  element  mesh  with  9 


106  M’.-//.  Clwn  Cl  til. 


KlilMFN’ 

D.O  F. 

3  0 

64 

1—211—1 

1 


I  1,1.  J  I  mile  cIcnu’Mi  mesh  Tor  ;i  ccniei-cmcket!  phitc. 


•/w 

/•■«/.  4.  V;iri:ilii>n  of  norni;ih/al  svinfiicirie  k-mlin^  Mress  inicnsily  f;it‘lors  for  cnick  Icnplhs  (nxuinplo  \i 

clcmcnls  arc  alsi)  shown.  'Ihc  varialion  of  synimciiic  bentiing  slrcss  intensity  factors  Kj 
normalized  by  the  value  of  hA/^  u/Zr  versus  ilie  crack  Icngiii  is  displayed  in  Fig.  4.  The  results 
arc  compared  with  the  solution  obtained  by  AVilson  and  Thompson  [8].  Good  agreement  is 
observed.  To  evaluate  the  accuracy  of  the  liniic  element  model  developed,  meshes  with  both  9 
and  25  elements  were  tested  and  very  little  difference  was  found.  Mcncc,  without  using  any 
singular  elements  or  .special  (inilc  element  models,  a  simple  mesh  with  9  elements  is  sumciently 
accurate  for  this  ease  even  as  the  ratio  of  w,  M’  reaches  O.S. 

The  rapid  convenience  of  the  bending  stress  intensity  factor  computed  can  be  seen  in  Table  I. 
Usually,  only  two  or  three  iterations  are  needed  to  obtain  acce|>tablc  results.  Table  2  presents 
the  comparison  of  the  polyiunnial  onler  N  taken  for  lining  Ihc  bending  moment  on  crack 
surfaces,  it  is  concluded  that  N—  I  or  2  is  sullicicnt  for  the.se  case.s. 
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Example  2:  A  rccuimpilar  cenier-cnu  kc.l  ihin  pU.ie  suhiecic.l  to  umforn,  pressure  will, 
various  Iwtindary  conditions 

The  geometry  ami  finite  element  mesh  are  similar  to  those  of  Rsample  I.  except  that  the  aspect 
ratio  of  the  plate  ll/W  is  taken  as  1.2.  Three  din'erent  houiulary  eomlitions.  i.c..  all  edges  simply 
supported  (SSSS),  two  edges  simply  supported  anil  the  other  two  edges  clamped  (SCSC).  and  all 
edges  clamped  (CCCC)  arc  solved  respectively.  The  variation  of  computed  symmetric  bending 
Stress  intensity  factors  versus  crack  lengths  for  the  SSSS  case  is  shown  in  Fig.  5.  For  comparison 
purposes,  the  analytical  solution  derived  by  Kccr  and  Sve  [261  and  hybrid  displacement  finite 
element  solution  computed  by  Chen  and  Chen  [  1 3]  arc  also  shown. 


Ttihh  I-  C*)nvcri!cncc  of  norni;ili/L*tl  iKmiing  stress  intcnxily  fatior  Ki 
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Fig.S.  Vari;ihoi»  o{  normalizcil  ^ylnmclr^c  iKiuhnji  sircss  iiilcnMly  facuuN  crack  icnplliN  umlcr  various  houiulary 

conditions  | Example  2|. 
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lixampiv  3:  A  nraimfuitir  thin  phiir  with  two  ntllincor  cracks  stthicclcil  to 
wtifonti  hcntlimi 

To  examine  ihc  imcraclion  cOccl  aition^  cracks  aiul  geometrical  Innmdarics.  the  problem  of  a 
thill  rectangular  plate  containing  two  collincar  inicmal  cracks  siibjecteci  to  unirorni  bending  at 
+  II  is  shown  in  ITg.  6.  I  he  dimensions  arc:  //  =  3\\\  11  =  lOii  and  h^a.  Kigure  7  displays 
variation  of  symmetric  Ixmding  stress  intensity  factors  normali/.ed  by  the  value  of 
versus  the  ratio  of  cfa%  where  c  is  the  half  ilistance  between  two  cracks.  A  linite  element  mesh 
with  9  elements  was  again  adoptetl  for  this  ease.  Ilowever.  as  cracks  approach  the  boundaries 
(e.g.  <■/</>  7).  the  residual  stress  near  the  boiindaiies  becomes  nu>re  complex  and  more  elements 
need  to  be  included. 


I  in.  u.  l  iiiiic  nxe^h  tor  a  plalc  willi  iwo  collincar  cracks. 


0  0  to  fO  aO  4  0  0.0  0.0  TO  oo 

e/« 

Hr/.  7.  Variation  »»r  norniali/cti  syni metric  Ix'iitiint:  stress  intciisily  factors  versus  the  ratio  v/a  (i’.xampic  3). 
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As  is  known  ihc  stress  distribution  in  the  vicinity  of  the  crack  lip  and  the  tending  sircsK 
iuisity  factor  depend  on  llic  plate  thickness  and  llie  specific  plate  theory  [3  5].  From  the  work 
rfHariranfl  and  Sih  [4],  it  is  found  that  the  present  computed  symmetric  tending  strew 
^ity  factor  K,  based  on  the  classical  plate  theory  can  be  deduced  from  that  derived  by 
Idssni’  theory  after  dividing  by  a  factor  <A(  I  )•  which  is  a  function  of  r  and  /./«.  Phe  value  of  the 
^;^|)  used  here  is  0.744578313  for  v=0.3  and  /.M=  1  [4]  Thus  for  --P”’ 
k  results  based  on  Rcissner’s  theory  obtained  by  Boduroglu  and  l-.rdogan  [71  for  a  plate 
■hounded  in  i«-direclion  are  calibrated  and  shown  in  l-ig.  7. 

TTie  interaction  cITcct  between  cracks  can  be  seen  for  the  ratio  (•/<;<  1.0.  The  tending  stress 
■tensity  factor  for  crack  tip  /1(Km)  is  larger  than  that  for  crack  lip  As  I.O«Ai<7.5.  as 

«oold  te  expected,  the  plate  with  two  collinear  cracks  behaves  like  nothing  but  an  inhnile  plate 
ootaininc  one  crack.  That  is,  neither  Ihc  interaction  elTecl  between  cracks  nor  he  boundary 
idluence  is  noted.  When  the  value  c/o  is  greater  than  7.5,  both  values  of 
«niHlv  tK  is  sliuhtly  larger  than  K,.)  due  to  the  influence  of  the  boundary.  As  c/a  ruichcs  8.0. 
Se  tecom«  e^  -o  baling  two  collinear  edge  cracks.  The  present  finite  element 
.negating  method  can  also  be  modified  to  deal  with  such  problem  by  extending  a  fictitious 
otek  which  is  the  same  as  the  original  crack  to  Ihc  outside  of  the  plate.  Again,  as  sa:n  in  Fifr  . 
J^lculated  result  by  this  modified  procedure  is  also  in  good  agreement  with  the  calibrated 

data  by  Boduroglu  and  Erdogan  [7], 


k  Conclusions 

Itased  on  the  analytical  solutions  derived  for  a  thin  cracked  plate  subjected  to  symmetric 
iSary  tending  on  crack  surfaces,  an  ellicient  and  accurate  finite  e  ement  allcrnaimg 
procedure  has  been  successfully  developed  to  deal  with  plate  tending  problems  with 
ll^ple  cracks.  Using  only  a  simple  convcniional  liiiite  elenicnl  mesh,  cxceficnt  f‘ 

jiromctric  tending  stress  iiilensily  factors  can  be  conipnteil.  1  he  mieraciioii  elTeci  .imoiit  cr.icks 
S  infitcMicc  tif  the  geometric  bou.ularies  on  the  calculation  of  the  teiidtng  stress  iiiteiisity 

factor  can  be  evaluated  accurately.  -  .  i  .1 

By  the  techiiiMue  developed,  the  computation  can  also  be  l^Tfonned  to  aiialy/e  the 
more  cracks  (c.g.  more  than  two  cracks)  witliout  any  ililliculty.  this  work  can  be  furlhtr 
atended  to  deal  with  thick  or  iiioilerate  thick  crackeil  plate  uiiilcr  arbitriiry  bending  type 
loadings.  This  will  be  presented  in  a  subsei|ueni  report. 
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